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BACKLASH  IN  SECON1  ORDER 
FEEDBACK  CONTROL  SYSTEMS 

1 .  Introduction 

Whenever  components  are  mechanically  connected  by  a  linkage 
or  a  gear  train,  there  exists  the  basic  problem  of  providing  smooth 
operation  without  looseness  in  the  connecting  parts  of  the  mechanical 
system.  In  general  this  basic  problem  is  not  completely  solved,  and 
the  connecting  parts  are  free  to  move  independently  over  a  limited 
range.  This  phenomenon  is  commonly  called  backlash  and  its  effects 
are  usually „  but  not  always,  undesirable.  In  particular,  when  back- 
lash  occurs  in  a  feedback  control  system  it  often  causes  a  small 
amplitude  "dither11  or  limit  cycle,  and  when  the  system  reaches  static 
steady  state  without  a  limit  cycle  the  accuracy  of  the  system  is 
affected  by  the  backlash. 

In  the  usual  feedback  control  system,  the  backlash  occurs  in 
two  places:  in  the  power  drive  gear  train  between  motor  and  load, 
and  in  the  gear  train  between  load  and  output  measuring  device.  It 
is  not  difficult  to  establish  the  pertinent  differential  equations  for 
either  of  these  cases,  nor  is  there  any  conceptual  difficulty  in  the  solu¬ 
tion  of  these  equations,  but  considerable  labor  is  required  for  the 
solution  of  even  the  simplest  cases.  Furthermore,  the  solution  of 
one  or  two  specific  cases  does  not  provide  perspective  for  generaliza¬ 
tions  s  and  the  multiplicity  of  parameters  involved  makes  detailed 
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analyses  unattractive.  This  report  deals  with  some  initial  results  in 
a  long  range  project  to  study  the  backlash  phenomena  quantitatively, 
with  the  ultimate  goal  of  providing  charts  which  will  be  useful  in  the 
engineering  analysis  and  design  of  feedback  control  systems  involving 
backlash  o 

In  general,  the  simplest  linear  feedback  control  system  has  a 
second  order  differential  equation,  and  is  commonly  referred  to  as  a 
second  order  system.  Backlash  may  and  does  occur  in  such  systems. 
While  it  is  obvious  that  most  feedback  control  systems  are  third  order 
or  higher,  the  study  of  backlash  in  this  report  is  restricted  to  second 
order  systems  only,  so  as  to  provide  a  firm  foundation  for  future 
studies  while  limiting  the  variable  parameters  to  a  variety  which  makes 
numerical  calculations  feasible.  In  like  manner  the  mechanical  para¬ 
meters  are  restricted  to  inertia  (mass)  and  viscous  friction;  couloumb 
friction,  though  of  some  practical  importance,  is  not  considered  quanti¬ 
tatively  and  is  permitted  in  a  qualitative  sense  in  only  one  special  case. 
It  is  recognized  that  the  mechanical  inertia  and  viscous  friction  may 
be  subdivided  and  appear  on  both  sides  of  the  backlash  in  varying  pro¬ 
portions,  such  cases  are  considered  in  detail.  It  is  also  recognized 
that  the  gear  tooth  surfaces  must  hit  each  other  after  the  meshing  gears 
have  gone  through  the  dead  space  (backlash),  and  this  impact  may  be 
either  plastic  or  elastic;  this  effect  is  also  studied  in  detail.  The 
mechanical  load  on  the  gear  train  is  considered  to  be  only  inertia  and 
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viscous  friction,,  This  combination  contains  only  one  energy  storage 
device  (inertia)  and  is  therefore  a  first  order  load*  In  some  cases 
the  mechanical  connection,  e.g.  shaft,  may  act  as  a  spring,  thus  pro¬ 
viding  a  second  order  load  which  is  capable  of  resonance  if  excited  by 
periodic  motion  of  the  proper  frequency.  This  case  is  not  considered 
but  is  reserved  for  future  study. 

The  two  basic  tools  for  studying  the  nonlinear  effects  of  backlash 
are  the  phase  plane-phase  space  technique  and  the  describing  function 
method.  The  describing  function  method  has  been  applied  to  a  number 
of  backlash  studies,  and  it  is  intended  to  devote  some  future  efforts  to 
an  extension  of  such  studies;  however,  the  investigations  reported  on 
here  are  based  solely  on  the  phase  plane-phase  space  technique.  For 
second  order  systems  the  phase  space  is  two  dimensional,  that  is, 
it  reduces  to  a  phase  plane,  and  nonlinear  as  well  as  discontinuous 
systems  can  be  analyzed  using  the  well  1  n  graphical  methods  of 
integration  (method  of  isoclines,  Lienard!s  construction,  the  Delta 
method)  which  make  the  phase  plane  a  powerful  tool,,  This  provides 
a  convenient  starting  point  for  the  analysis,  though  the  graphical 
methods  do  not  apply  to  higher  order  systems. 

The  phase  plane  technique  with  graphical  integration  was  used  in 
initial  studies  recorded  here  but  was  applied  to  only  the  simph 
cases.  The  results  obtained  and  the  procedures  developed  provided 
sufficient  insight  to  permit  formulation  of  a  rather  comprehensive 
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digital  computer  program  to  study  the  effects  of  many  parameters  on 
the  response  of  a  given  system.  The  majority  of  the  data  presented 
here  was  obtained  using  a  digital  computer  program;  however,  any  of 
the  results  so  obtained  could  be  duplicated  by  graphical  methods. 

The  digital  computer  was  programmed  to  provide  a  complete  solu¬ 
tion  of  the  differential  equations,  and  sufficient  data  was  computed  to 
provide  a  complete  and  very  accurate  phase  trajectory.  However,  for 
most  of  this  initial  phase  the  question  of  primary  interest  was  the 
development  of  criteria  for  the  existence  of  limit  cycles,  The  computer 
program  was  therefore  provided  with  a  subroutine  which  printed  out 
only  the  maximum  position  per  cycle.  This  was  done  to  decrease  the 
read  out  and  print  out  times.  By  a  simple  auxiliary  command,  complete 
data  for  the  phase  trajectory  was  made  available;  however,  the  com¬ 
plete  phase  trajectory  was  solved  for  only  a  restricted  number  of  cases 
because  of  the  time  involved.  As  a  result,  the  limit  cycle  studies  are 
detailed  and  fairly  complete  (more  data  would  improve  the  results  but 
should  not  introduce  any  startling  features);  but  many  questions  con¬ 
cerning  response  during  the  transient  period  have  not  been  completely 
answered.  Much  of  the  remaining  work  (within  the  restricted  scope 
of  this  initial  phase)  is  just  a  matter  of  data  reduction. 
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2 .  THE  NO  LOAD  CASE 


For  some  systems  one  may  neglect  inertia  and  viscous  friction 
on  the  output  side  of  the  gear  train  containing  the  backlash.  This  is 
arbitrarily  called  the  No  Load  Case.  Two  practical  variations  are 
common.  First,  the  most  commonly  encountered,  is  the  case  of  the 
instrument  servo  which  is  driving  some  simple  indicating  or  measur¬ 
ing  device  such  as  the  wiper  arm  on  a  potentiometer.  In  this  case 
the  actual  driven  load  has  a  negligibly  small  inertia  and  viscous  fric¬ 
tion.  There  is  usually  some  couloumb  friction  in  this  load,  but  it  is 
often  too  small  to  affect  the  differential  equations  though  large  enough 
to  prevent  jitter  in  the  low  inertia  driven  load.  The  second  case  of 
interest  here  is  that  in  some  larger  systems  with  driven  load  of  appre¬ 
ciable  inertia 5  the  power  gear  train  has  negligible  backlash,  but  there 
is  a  measuring  device  geared  to  the  load,  and  there  is  backlash  in  the 
gear  train  between  load  and  measuring  device.  Note  that  while  these 
two  cases  arise  in  systems  which  are  physically  different,  they  are 
mathematically  equivalent  and  thus  one  analysis  suffices. 

The  block  diagram  description  of  the  system  is  given  in  Fig,  2-1. 
11  it  is  assumed  that  a  step  displacement  input  is  applied  to  the  system 
when  the  gears  are  in  contact  and  the  system  is  driving  linearly,  then 
the  basic  linear  differential  equations  are. 


E=eR-  ®C 
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J 


A 

Backlash 


Fig.  2-1  Block  Diagram 


r 


The  dynamic  response  of  the  system  is  shown  on  the  phase  portrait 

of  Fig.  2-2,  and  equations  2-1  and  2-2  describe  the  first  portion  of  the 

trajectory  from  M  to  N.  At  point  N,  E  =  O  and  ©  =  0?  the  motor  re- 

verses  and  runs  in  the  opposite  direction,  increasing  the  velocity  and 

displacement  of  the  motor  shaft.  However  the  output  shaft  cannot  move 

and  the  error,  E.,  cannot  change  until  the  backlash  is  taken  up.  The 

operation  during  this  interval  is: 

JM0M+fM@M  =  KEN  (2-3) 


from  which 


where 


•f=f/2VKJ  M  andW" 


(  2-4) 

(-  l+2fWrt+e'2^Wi'1r  ) 


Also 

’  _  -2fWr»t 

0M--"2^(1-e  ) 


(2-5) 


It  is  readily  seen  that  the  horizontal  or  E-axis  is  a  locus  of  points 
at  which  the  describing  differential  equation  changes  discontinuously , 
and  it  may  then  be  called  a  "dividing  line."  This  particular  dividing 
line  is  a  "separation  line"  i.e.  as  the  state  point  crosses  this  line  the 
physical  system  goes  into  the  backlash  region  and  the  drive  member 
is  mechanically  separated  from  the  load  member  During  the  time 
interval  when  equations  2-4  and  2-5  apply,  and  the  backlash  is  being 
"taken-up";  the  load  member  does  not  move,  E  and  E  do  not  change, 
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FIG.  2-2  PHASE  TRAJECTORIES  OF  AN  INSTRUMENT  SERVO¬ 
MECHANISM  WITH  BACKLASH 
Di  -DIVIDING  LINE  FOR  RECOMBINATION 
D2  -DIVIDING  LINE  FOR  SEPARATION 


but  ©  and  ©  do  change:  On  the  phase  plane  of  Fig  2-2  the  vertical 
mm  ° 


line  from  N  to  P  indicates  the  first  interval  of  this  type.  It  can  be 
argued  whether  this  line  is  or  is  not  validly  a  segment  of  the  phase  tra¬ 
jectory,  but  the  physical  interpretation  is  quite  clear:  E  and  E  are 
constant  at  point  N  until  the  backlash  is  completely  taken  up,  then  as 
the  gear  teeth  hit  (with  inelastic  impact)  the  load  instanteously  acquires 
a  velocity  and  the  state  point  jumps  discontmuously  from  N  to  P  Once 
the  backlash  is  taken  up  the  system  once  more  operates  linearly  and 
the  segment  of  phase  trajectory  from  P  to  Q  is  again  governed  by 
equations  2-1  and  2-2.  The  point  P  is  one  point  on  another  dividing 
line,  this  is  the  line  at  which  the  backlash  is  taken  up  and  the  system 
physically  recombines,  i.e.  ,  the  mechanical  connection  is  restored. 
Thus  the  line  may  be  called  a  backlash  n recombination  line0n  This 
line  can  be  constructed  accurately  only  by  a  point-by-point  calculation. 
However,  equation  2-5  can  be  rewritten: 


(2-6) 


N 


where  is  the  time  interval  required  to  take  up  the  backlash  For 


values  of  the  backlash,./!^  ,  the  recombination  line  may  be  represented 
on  the  phase  plane  as  shown  in  Fig.  2-3.  Only  the  tangent  line  is 
accurate  (unless  detailed  computations  are  used),  but  knowledge  of 
this  tangent  line  permits  derivation  of  a  stability  criterion. 
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FIG.  2-3  FAMILY  OF  BACKLASH  DIVIDING  LINES  WITH 
COMMON  TANGENT  AT  THE  ORIGIN. 

.4 MAX*  TAN'1—  1/2 j 
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2 . 1  A  Stability  Criterion  for  the  No  Load  Case 


If  equations  2-1  and  2-2  are  combined  and  nondimensionalized 
forming  a  single  equation  with  error,  E,  as  the  independent  variable, 
this  equation  can  be  manipulated  and  integrated  to  define  the  linear 
phase  trajectory  on  the  phase  plane  as* 


a  2  #  2  2 
E  +2f  wnEE  4*  UJ  f\  E 


2iT 


-1 


tan 


2  V  1— ? 

=  D  e 


E+f  O  n  E 
WnE-v/  1-f  * 


Using  the  linear  transformation 


(2-7) 


Y=E+£wnE,X=  wnE  —  £ 


2  2  2 
X  +Y  =D  e  v 


2  JT  .  -1  7 

- tan  — 


(2-8) 


(2-9) 


-1 


Letting  Y=  j>  sin  yp  ,  X=  ^  cos  yp  ,  \p  =  tan  Y/X 

Tv* 


(2-10) 


p  =  De 

The  backlash  dividing  lines  are  located  on  the  X  -  Y  plane  by  noting 


that 


Y 

X 


O 

_  E+jTa»n  E 
UpE^l-  f  2 


(2-11) 


but  the  dividing  line  for  separation  is  the  E-axis,  for  which  E=Os  and 
for  this  line 


Y_ 

V 
u.  V 


r 


v 


i 


J* 


2 


tan  yj, 


(2-12) 
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The  backlash  dividing  line  for  recombination  is  not  readily  determined, 


but  the  tangent  to  this  curve  A  on  Fig.  2-3,  is  easily  transferred: 

max 


Y_  E/  uy 
X  E 


(2-13) 


s/i-r2  V i-  $ 


and  for  E  /  E  =  —  1  /  2 

max  ° 


so 


Y  2^-1 


X 


2fV  1— 


=  tan 


(2-14) 

max 


which  may  be  expressed  as  an  acute  angle 


tan 


max 


1  —  2f 2 

2|\/l_  i-2 


(2-15) 


In  Fig.  2-4  the  backlash  line  for  recombination  is  sketched,  not  calcu¬ 
lated*  but  the  basic  condition  for  increasing  oscillations  is 

Pl"1  'P  1  =  P2  cos  N't  ,2-16) 

Assuming  that  the  angles  are  measured  from  the  negative  real  axis  in 


Fig  2-4  and  that  p 

and  are  known: 

£  $  1 

D=pie  \/  1  —  i*2 

(2-17) 

J*(*P2-  1  ) 

„  Vi—  J2 

Vi—  ?2 

(2-18) 

p2=De  =  pi* 


N>2=  7t  -  N>t 


(2-19) 
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Substituting , 


(2-20) 


cosl|^=-  cos 


V,.f2 


V 


/ 1 


(2-21) 


The  angle  may  be  chosen  arbitrarily,  and  by  inspection  of  Fig.  2-4 
it  is  apparent  that  if  the  required  inequality  exists  it  must  be  most  pro¬ 
nounced  when  the  point  (p^  ,  1)^)  is  l°cated  near  the  origin  on  the  tangent 

line,  in  which  case  =  U>  Substituting  this  value  for  UJ  * 

1  1  max  1  1 


(2-22) 


l£—  e 


f  - 1  i  -  2  i*2  ^  -1  £ 

•  —  tan  — - -  -  1 —  —  tan  - 


Vi-/ 


2fVl-  t 


\/l  —  t 


Equation  2-22  has  been  plotted  in  Fig.  2-5,  and  the  stability  point  is 
marked.  This  shows  that  the  system  with  backlash,  but  no  load,  will  have 
a  limit  cycle  if  £  0.29,  and  the  existence  of  the  limit  cycle  is  inde¬ 

pendent  of  the  amount  of  backlash,  so  long  as  some  backlash  is  present. 
Obviously  the  amplitude  of  the  limit  cycle  will  be  determined  by  the  amount 
of  backlash. 
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X|'P|  COS  ¥|  ,  X2  =  02C0S'i'T 

*T  *  tan-2  c/yi-T2 

T— TANGENT  TO  DIVIDING  LINE  AT  ORIGIN. 
L— SEGMENT  OF  LOGARITHMIC  SPIRAL. 


0  .1  .2  .3  .4  .5  .6 

FIG.  2-5  GRAPH  OF  STABILITY  CRITERION  FOR  BACKLASH. 

k  =  e  ~a  (n- tan*'  b-  tarf1  a ) 

2C 

a'c/yr-s2 

b=(l-252)/2jy]-52  15 


3.  PHASE  PLANE  ANALYSIS  C  F  BACKLASH  IN  A  SECOND  ORDER 
SYSTEM  WITH  DIVIDED  LOAD  AND  PLASTIC  IMPACT. 

In  the  more  general  cases  of  backlash,  the  inertia  and  viscous  fric¬ 
tion  on  the  output  side  of  the  gear  train  are  not  negligible  .  For  second 
order  systems,  such  cases  can  be  analyzed  using  phase  plane  techniques. 
This  is  true  for  any  subdivision  of  inertia  and  friction  between  drive 
member  and  driven  member,  for  any  gear  ratio,  and  for  plastic  or 
elastic  impact  at  the  gear  tooth  surface.  While  actual  graphical  integra¬ 
tion  to  construct  the  phase  trajectories  is  too  laborious  to  be  a  tool 
for  generalized  analysis,  understanding  the  techniques,  and  actual  analysis 
of  one  or  two  specific  cases  helps  build  appreciation  for  the  nature  of  the 
problem  Graphical  solution  points  out  the  physical  phenomena  and 
boundary  conditions,  and  leads  naturally  into  the  preparation  of  digital 
computer  programs. 

When  the  backlash  is  taken  up  and  the  system  is  operating  linearly, 
the  equations  of  motion  are  exactly  the  same  as  for  the  no  load  case. 
Considering  here  that  the  motor  used  is  specifically  a  D.  C<  armature 
driven  motor,  and  using  a  detailed  notation,  the  equations  of  motion  are: 


(  n2  j  +  jt  )  e  +—  (N2  f  +  f  )  ©  =  t 

N  '  m  L  c  N  m  L,  c 


T=  K1  I 

o 

V— R  I  +  K  e 

2  m 

V  K2 

i  =  (TT  —  N  ec> 


(3-1) 

( 3=2) 
(3-3) 

(3-4) 
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(3-5) 


V  =  K  (  ©  -  ©  ) 
3  K  c 


4"  (  N2  J  +  JT  )  ©  +4?  (N2  £  +  f_  + 

N  m  L  c  N  m  L 


•S  v-2 


N 


R 


)  © 


K  K  © 

+  1  3  c 

R 

Where 


K1  K3  ®R 

R 


( 3-6) 


K=  Motor  torque  constant 

1=  Motor  generated  voltage  constant 
Error  measurement  constant 
■=Kj  K3/R 
R=armature  resistance 
N=gear  ratio 

J=  inertia  of  system  without  load 
J1?  load  inertia 

£=  viscous  friction  (mechanical)  of  system  without  load 
f=  viscous  friction  of  load 

F=  N  f  +  K_  K  N/R  =  viscous  damping  of  system  without 
m  m  12/  ,  ,  ® 

load  but  including  motor  back  EMF 

5:£l/n 

^commanded  displacement 
<3=  displacement  of  motor 
(^displacement  of  load 
©^displacement  of  combined  system. 

To  permit  convenient  numerical  manipulation,  let  N=  1  and  0  =  1, 

R 

then  equation  3-6  becomes: 


•  *  p*  p* 

©  +  (  m  +  L  )  ©  =  ( 

C  c 

J  +  JT 
m  L 


K  )  (1  -  ©  ) 

c 


(3-7) 


J  +  JT 
m  L 


Equation  3-7  may  be  manipulated  into  the  isocline  equation  in  the  usual 


fashion: 


©  =  1  -  © 
c  c 


N  ,+ 


F  Ft 
m  +  L 


1  J  +  J 
m  L 


K/  (  J  +  JT) 
_  /  m  L 


(3-8) 


where  N  is  defined  to  be  d  ©  /d  ©  and  is  the  slope  of  the  phase 
1  c/  c 

trajectory  as  it  crosses  the  designated  isocline. 
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Equation  3-8  is  the  normal  isocline  equation  for  a  linear  second  order 


feedback  control  system.  It  is  expressed  in  terms  of  motor  inertia 
and  friction  and  the  load  inertia  and  friction  because  this  subdivision 
becomes  important  when  backlash  is  considered 

When  there  is  backlash  between  motor  and  load,  and  when  this  back¬ 
lash  is  inside*  the  closed  loop,  the  system  enters  the  backlash  region 
(i.e.  ,  the  gears  lose  contact)  when  the  slope  of  the  load  deceleration 
trajectory  is  exactly  equal  to  that  of  the  linear  system.  In  terms  of 
equations,  when  the  load  is  drifting  free- 


J 


L 


(3-9) 


from  which 


d© 


L 


(3-10) 


When  the  system  is  operating  linearly  as  described  by  equation  3-8,  and 

the  state  point  reaches  the  isocline  which  is  numerically  equal  to 

then  the  gears  lose  contact,  l  e.  ,  the  backlash  separation  line  is  that 

isocline  of  the  linear  system  for  which  N  =  N  =  —  F  I  J 

1  L  L  /  L 

When  the  backlash  separation  line  is  crossed  the  load  drifts  freely 
according  to  the  relationship  of  equation  3-10,  which  guarantees  a  phase 
trajectory  that  is  a  straight  line  The  motor  shaft  does  not  remain 
stationary  It  moves,  but  is  driven  open  loop  rather  than  closed  loop, 


*It  is  possible  to  have  backlash  in  the  gear  train  between  motor  and 
load  but  so  located  that  it  is  between  the  measurement  pick  off  and  the 
load.  Then  there  is  no  backlash  in  the  closed  loop,  which  operates 
linearly.  This  case  cannot  have  a  limit  cycle  and  is  not  considered  here. 
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with  a  time  varying  input.  The  equations  for  the  motor  shaft  motion  is 
•  «  « 

J  ©  +  F  ©  =K(©o  —  ©.)  (3-11) 

m  m  m  m  R  L 

The  phase  plane  representation  of  these  equations  is  shown  in  Fig,  3-1. 
The  separation  dividing  line  is  the  isocline  passing  through  21.  Equation 
3-11  cannot  be  represented  by  an  isocline  relationship,  so  first  the 
analysis  requires  a  time  solution  of  equations  3-9  and  3-11,  The  back¬ 
lash  is  completely  taken  up  in  the  reverse  direction  when  the  difference 
between  the  displacement  of  the  load  and  the  displacement  of  the  motor 
is  exactlvA,  the  magnitude  of  the  backlash.  In  the  solution  of  equations 
3-9  and  3-11  the  initial  conditions  are  determined  at  the  separation  dividing 
line  and  are  as  shown  on  Fig.  3-2.  Transforming  equation  3-9 

=  O  (3-12) 


2 

s  ©T  (s)  —  s  ©  — -  © 

+  F 

s  e  (s)  —  0, 

L  Li  h 

L 

_  ]} 

L 

from  which 
eT  (s)  = 


© 


h 


L'"'  s(s+Fl/Jl) 


(3-13) 


Taking  the  inverse  transform 


fi/ 


e  (t)  =  e  + 

L  L.  FLjJL 

Equation  3-11  transforms  to 


J 

©  (s)  —  s  ©  —  © 

+  F 

s  ©  ( s)  —  © 

m 

m  m  ^  m^ 

m 

m  m  ^ 

=  K 


(3-14) 


- ©  (s) 

S  L 


(3-15) 
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INERTIA  AND  FRICTION  ON  BOTH  SIDES  OF  BACKLASH 


N  =  -  f/J  =  SEPARATION  ISOCLINE;  A=  TRAJECTORY 
FOR  COMBINED  SYSTEM;  B= TRAJECTORY  FOR 
DRIFTING  LOAD;  C=  TRAJECTORY  FOR  MOTOR 
AFTER  SEPARATION. 
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I  T*  ©Ml  I  —  ^Ll 

TAN¥  =  s  - 6 


ec*eR«i.o 


IG.3-2  RELATIONSHIP  BETWEEN  9U  ,  6 


since  ©  (t)  =  u  (t)  and  ©  (s)  =  l/s. 

£\  K  ^ 


R 


Substituting  equations  3-13  for  ©  (s),  and  noting  equivalent  initial 

conditons  from  Fig.  3-2,  equation  3-15  may  be  manipulated  to 


m 


©  (s)  - 

m 


©  K{  1  -  ©  ) 

m  j  mi 

hs(  s  +F  |J  )  +  1  s^(s+  F  |J  ) 

mm7  m  ml  m 


K  © 


m . 


J  s  (s  +  F_  J_  )  (s  +  F  I  J  ) 
m  L  L  m 


)m 

Taking  the  inverse  transform  of  equation  3-1  6  gives: 

J  ©  Ft  K  (1  -  ©  )  J 

mm  — 

©  (t)  =  ©  +  — - -  (  1  -  e  J 

m  m  F 

1  m 


m 


m 


)- 


m  '  m 


"F  T 

m 


( 3-  1 6) 


(  1 


F  t 
m 
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F  F 
m  L 


-  e 


F  t 
m 


m 


J  F  -  J  F 
L  m  m  L 


*  ( 


m 


K  © 


m 


)  - 


F  F 
L  m 


F  t 
m 


m 


r-  J  J 
m  _  L 
t  _  —  —— 

m  L 


F  1? 


(3-17) 
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m  L  L 

If  equations  3-14  and  3-17  are  plotted  with  displacement  as  ordinate 


and  time  as  abscissa  they  may  be  used  as  a  nomograph  to  evaluate  points 
on  the  backlash  dividing  line  for  recontacting  of  the  gear  teeth.  If  both 
are  plotted  on  the  same  coordinate  axes,  then  at  any  instant  the  difference 


between  the  ordinate  values  is  the  relative  displacement  ©  -  ©  , 

Lj  m 

and  whenever  ©  -  ©  =  j£\  the  backlash  is  taken  up  Since  this  point 

L  m  ^ 

also  defines  values  for  ©  and  for  t,  differentiation  of  the  proper 

L 

equations  permits  evaluation  of  ©  and  ©  at  that  instant. 

m  1j 
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For  repeated  cycles  these  time  curves  need  not  be  replotted,  but  may 


be  used  as  a  nomograph  by  simply  adjusting  the  ordinate  scales.  This 
is  due  to  the  fact  that  the  only  change  in  the  physical  situation  is  a  change 
in  the  initial  conditions,  and  from  Fig.  3-2  there  is  a  fixed  relationship 
such  that  i  q  1  —  © 

tan  t|>  =  -  =  --g -  (3-18) 


ml 


1 


where  tan  4*  is  known  number.  Then  1  -©  =  ©  tan  4*  which  (3-19) 

m  j  m  j 

indicates  the  required  adjustment  in  ordinate  scales.  A  pair  of©  and 

Li 

©  curves  is  shown  on  Fig.  3-3. 


m 

At  the  instant  backlash  is  taken  up  both  motor  and  load  have  finite 
velocities  .  Therefore  both  the  law  of  conservation  of  momentum  and 
the  law  of  conservation  of  energy  must  be  considered  If  the  impact  is 
defined  to  be  inelastic  (plastic)  this  is  equivalent  to  stating  that  sufficient 
deformation  occurs  to  dissipate  some  energy  as  heat;  then  the  law  of 
conservation  of  energy  need  not  apply  to  the  dynamics  of  the  mechanical 
system,  but  the  law  of  conservation  of  momentum  still  restricts  the 
boundary  conditions.  It  follows  that  after  recontact  of  the  gear  teeth, 
the  requirement  that  momentum  must  be  conserved  is 


o  o 

j  ©T  +  j  © 

L  L  mm 


(  J  + 
m 


V 


(3-20) 


where  is  the  instantaneous  velocity  of  the  recombined  system.  Note 
that  the  points  at  which  the  gear  teeth  recontact  define  a  dividing  line  on 
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the  phase  plane,  and  the  result  of  the  applying  the  momentum  restoration 
defines  another  dividing  line  on  the  phase  plane.  These  lines  may  be 
established  point  by  point  in  evaluating  the  phase  trajectory,  but  no  simple 
equation  for  these  lines  has  been  obtained. 

A  phase  plane  plot  for  a  particular  case  is  given  on  Fig  3-4  This 
plot  shows  the  three  dividing  lines,  and  the  nomograph  of  Fig.  3-3  was 
used  in  computing  Fig.  3-4.  Note  that  the  system  goes  into  a  limit  cycle 
along  path  ABCD, 
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Combined  System  0c  + .0.2 0c  +  0c  =  0r 

System  with  load  Separated  0.5©m  +  0.2©m  =  0r-0l 
Load  Only  0.5  ©l  =  0 

Backlash  =  0.3  RAD 


4.  COMPUTATION  BY  DIGITAL  COMPUTER 


4.  1  Introduction 

Use  of  the  phase  plane  techniques  of  section  3  is  instructive  but 
prohibitively  laborious  if  a  large  number  of  cases  is  to  be  investigated. 
Alternatives  are  analog  simulation  or  digital  computation.  The  indivi¬ 
dual  equations  are  readily  simulated  on  the  analog  computer,  but  instru¬ 
mentation  of  the  boundary  conditions  can  be  difficult.  An  analysis  using 
analog  simulation  has  been  carried  out  elsewhere^  but  was  restricted 
to  the  inelastic  impact  case.  It  was  decided  to  use  a  digital  computer 
to  calculate  the  desired  data,  and  the  computer  was  to  be  programmed 
using  phase  plane  concepts  so  that  all  of  the  data  required  for  phase 
plane  plots  would  be  available  when  needed.  The  first  program*  was 
prepared  for  the  inelastic  impact  case  only.  It  operated  successfully 
as  verified  by  point  checks  calculated  long  hand  on  the  phase  plane  . 

While  the  results  were  accurate,  only  a  limited  amount  of  data  was 
computed  due  to  lack  of  available  time .  Upon  resumption  of  the  pro¬ 
ject  it  was  decided  to  reprogram  the  computer  to  handle  the  general 
case  of  elastic  *  contact  of  the  gear  teeth,  including  the  case  of 
plastic  impact  as  a  limiting  case  This  program  was  completed 
and  operated  satisfactorily.  It  is  given  in  detail  in  the  appendices. 

The  results  reported  in  the  remainder  of  this  text  is  based  on  com¬ 
putations  made  with  this  program. 

♦Knoll  &  Narud 
*New 

^Anderson  &  Luckett 
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TABLE  4=1  SYMBOLS* 


®R 

(ONE)** 

Angular  step  input  (radians). 

Desired  load  displacement  at  steady 
state 

© 

c 

(THETAL) 

Output  displacement  of  load  (radians). 

0 

© 

c 

( THE  TAD  L) 

Output  velocity  of  load  (radians) 

sec 

E 

Error  in  radians.  E=©_  —  © 

R  c 

K 

(KMOT  CONST) 

Constant  of  proportionality  relating 
position  error  to  torque  developed  by 
the  motor . 

T 

m 

Developed  torque  of  motor. 

tl 

Torque  of  load. 

TLm 

Torque  of  load  referred  to  motor. 

P 

(RHO) 

Number  of  teeth  on  gear  1 

Number  of  teeth  on  gear  2 

J 

m 

(MOINERT) 

Inertia  of  motor  measured  at  motor 

]L 

( J  LOADVAL) 

Inertia  of  load  measured  at  load. 

J 

Total  inertia  measured  at  motor 

J  =  Jm+p^J 

L 

F 

m 

(FMOTOR) 

Friction  of  motor  measured  at  motor. 

fl 

(FLOAD) 

Friction  of  load  measured  at  load. 

A 

(DELTA) 

Backlash  of  gearbox  measured  in  radians 
at  the  output 

N 

s 

Slope  of  phase  trajectory  (combined 
system. 

d©  © 

XT  A  c  c 

NS  =  d©"  =  S~ 
c  c 

*In  this 

table  some  of  the  symbols  are  redefined  and  have  different 

meanings  than  in  preceding  sections. 

**(  )  Te 

rms  indicate  computer 

mnemonics  referred  to  in  Sec.  3, 

Computer  Program  Development. 
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nl 

Slope  of  phase  trajectory  (load  floating 
free). 

d©  © 

N  =  c  .  c 

L  d©  © 

c  c 

.1  ,1 

©  ,  © 

Velocities  of  motor  and  load  after  impact 

me 

e 

(RESTITUT) 

Coefficient  of  restitution. 

r 

(ZETA) 

System  damping  coefficient. 

(OMEGAN) 

System  natural  frequency. 

Wn 

(OMEGANSQ) 

System  natural  frequency,  squared. 

(FLFTPRIN) 

TT 

Load  riction  ratio  with  respect 

*  T 

to  friction  of  the  system. 

J 

m 

(INERTRAT) 

Inertia  ratio. 

JL 

© 

m 

(THETAM) 

Displacement  of  motor  (radians)  . 

© 

m 

(THETADM) 

Velocity  of  motor  {radians] 

V  sec  / 
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4  c  2  Development  of  Equations 


The  equations  developed  were  those  of  a  phase  plane  analysis. 

Prior  to  defining  the  net  system  equations ,  the  assumptions  under  which 
the  analysis  was  made  will  be  stated.  When  necessary  these  assump¬ 
tions  will  be  amplified  and  referred  to  later  in  this  work 
It  was  assumed  that: 

1  .  The  gear  teeth  were  initially  in  contact  and  the  initial 

conditions  of  the  system  were  all  equal  to  zero.  This  was 
later  proved  to  be  an  unnecessary  limitation  for  the  study 
of  steady  state  response. 

2.  Plastic  deformation  of  the  gear  teeth  during  steady  contact 
and  impact  and  any  torsional  deformations  of  driving  shafts 
are  negligible . 

3.  The  inertias  of  the  gears  and  drive  shafts  are  considered 
as  part  of  the  load  or  motor  inertia  depending  on  their 
attachment  in  the  system. 

4.  The  law  of  conservation  of  energy  was  completely  satisfied 
in  only  the  perfect  elastic  case  by  maintaining  the  total 
mechanical  rotational  energy  of  the  system  constant  at  the 
instant  prior  to  and  after  impact.  The  law  of  conservation 

of  momentum  is  satisfied  in  plastic,  elastic,  and  intermediate 
cases.  When  the  law  of  conservation  of  momentum  is  the  only 
equation  required  to  be  satisfied,  the  energy  lost  from  the 
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system  is  dissipated  in  the  heat  of  infinitesimal  deformations 


of  the  gear  teeth 

The  gear  teeth  are  in  contact  only  instantaneously  during 
impact  for  the  elastic  and  intermediate  cases  (excluding 
plastic)  and  that  the  impulse  torques  of  drive,  friction  and 
bearing  supports,  etc.  ,  are  zero  during  impact,  e.g. 

5 


T  dt  =  O  and  ©  =  0 
m 


dt  o 


6.  The  coefficient  of  restitution  of  the  two  opposing  gear 

teeth  is  the  same  or  is  described  by  an  equivalent  coefficient 
if  the  two  gear  teeth  are  of  unequal  coefficients. 

7.  Backlash  is  assumed  to  be  equal  at  all  points  on  the  gear 
circumference.  Backlash  is  measured  at  the  output  shaft 

A  block  diagram  of  the  system  considered  is  presented  in  Fig.  4-1, 
and  symbols  have  been  defined  in  Table  4-1. 

The  equations  for  the  motor  and  load  in  contact  are* 


<D 

1 

i 

® 

n 

H 

(4-1) 

T  = 

R  c 

KE  =  K  (eD  — 

©  ) 

(4-2) 

m 

R 

c 

T  = 

O  G  O 

JT  ©  +  fT  © 

(4-3) 

L 

L  c  L  c 

TLm 

=  P  Tl=P  <JLec+  fLec' 

(4-4) 

©  = 

O  ©  forward 

©  =  o  ©  +  A  backward 

(4-5) 

c 

r  m 

drive 

c  t  m  drive 
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Fig.  4-1  Block  diagram  of  second  order  servo¬ 
mechanism  with  backlash  in  the  gears  between 
motor  and  load. 


32 


0  =  p  0  forward  and  backward  drive 

c  r  m 

•  «  o 

hence  T  =  p  (J  ©  +  f_  p  ©  ) 

Lm  r  Lp  m  L  “  m 

=p2  (J.  ©  +  f  ©  ) 
r  L  m  -L  m 


T 

=  J 

g  o 

© 

0 

+  f  © 

+  Tt 

m 

m 

m 

m  m 

Lm 

J 

*  0 

© 

f 

2 

0  p 

(JT  ©  f.  © 
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c 

j-m 

L,  c  L  c 

p 


(j 


m 


p  j.)  ©  (f 
r  L  c  +  J2  + 
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P  9 

P  fL)  Gc 


T  =  K  (©  —  ©  ) 

m  R  c 


K©  =( - — +pJT)  ©  +  (—  +  pfT  )  ©  +  K© 

R'p  r  L  c  p  r  L  c  c 
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©  + 


(g  tPV-  e  + 


K  © 


C  (Jm  PJL) 

- + 
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k©t 


Jm  P J: 


(4-6) 

(4-7) 


(4-8) 


(4-9) 

(4-10 


(4-11) 


<•  h.  ‘  2  2 

©  +2YUL©  +  u>n  ©  =u)n  ©T 

c  V  0  c  1  *  C  1 


where  u>n  =  pK 


jm+p2j 


2  £<*>„= 


m  + 


a 


j  +  p  j. 


m 


(4-12) 

(4-13) 


(4-14) 
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The  equation  for  load  alone  with  gears  not  in  contact  is' 


J  ©  +  fT  ©  =  O 

L  c  L  c 


(4-15) 


or 


©  + 


c  J . 


©  =  O 

c 


The  equation  for  the  motor  alone  with  the  gears  not  in  contact  is. 

J  ©  +  f  ©  =  K  (  ©_  —  ©  )  (4-16) 

mm  mm  R  c 

The  method  of  analysis  used  by  New,  Ref.  e,  is  used  in<  order 
to  be  able  to  examine  the  system  independently  of  its  Ulp 
Defining  first  U)^t  =  t# 


and  differentiating  Lc/r)dt  =  dt* 


.  ,  dt*  ,  ,  ,  2  .  dt*  . 

U)  n  =  — —  and  =  (  — — — ) 


dt 

d© 


©  = 


d©  dt* 
c 


dt 
d© 


c  dt 


dt*  dt  ’  dt* 


c  ** 
=  © 


* 


hence  ©  =©  UJ 


n 


similarly 


d2© 


©  = 

c 


dt 


d2© 


dt  *1 


dt 


2 

d  © 


UV 


d2  © 


(at*) 

finally 


•  •?]>  «•  *  2 

-  =  ©  and  ©  =  ©  UJ  i 

2  c  c  c 


'n 


♦  •*  dt* 

9  =  ©  U>n=  ©  ~rr — 

c  c  c  dt 


Jecd‘=j 


©  dt  or  ©  =  © 

c  c  c 


(4-17) 

(4-18) 

(4-19) 

(4-20) 

(4-21) 

(4-22) 
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Using  the  equations  of  the  combined  system  and  the  load  alone, 

w  0  2  2 

©c+  2  y 0Jn©c  +  ©^  =  vjU n  ®r  combined  system  (4-23) 

£l  - 

©  +  — T  ©  =  O  Load  alone  (4-24) 

c  u  c 

L 


and  making  the  indicated  substitution  to  a  transformed  (*)  coordinate 
system  for  the  system  equation 

UJn  ec+2fu;n  ©c+u;n©c=Uin  ©R  (4-25) 

®c  +  2-^e  c  +  ®c  =eR  <4'26> 

and  introducting  the  slope  equations  of  the  phase  plane 


N  ©  +  2  {*  ©  +  ©  =  ©T 

g  n  +*  n  r-  T 


c  c  c  R 

If  the  load  equation  is  first  put  in  the  form  for  the  phase  plane 


(4-27) 


N  ©  +  qL  ©  =  O 

L  c  J  c 
J— i 

and  then  transformed  to  the  (#)  coordiante  system  where 


(4-28) 
©  £  O 


L 

N  ©  +  — —  ©  =  O  (4-29) 

Li  C  J  C 

L 

The  same  result  may  be  obtained  by  setting  ^n=  1.  The  results  of 
such  a  transformation  require  inverse  scaling  for  practical  applica¬ 
tion,  examples  of  which  given  later.  It  is  noted  that  the  slope  of  the 
load-free  equation  is  mathematically  the  same  whether  in  the  phase 
plane  or  in  the  transformed  phase  plane.  It  is  pointed  out  at  this 
time  that  the  equations  later  developed  to  satisfy  the  laws  of  con¬ 
servation  of  momentum  and  energy  are  independent  of  the  system 
natural  frequency. 
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To  establi  sh  equations  for  the  law  of  the  conservation  of  momentum 

and  a  relationship  satisfying  the  law  of  the  conservation  of  energy, 

the  schematic  of  Fig  4-2  is  used.  It  is  seen  that 

*  •  1 

©c=  —  p  © j  >  ©^-  —  p  ©i  and  since  p  =  1  an  equivalent 

schematic  is  shown  in  Fig.  4-2b.  By  representing  the  inertia  of  the 
load  and  motor  as  inertia  of  the  gears  figure  4-2c  is  obtained  at  the 
instant  of  impact.  The  torque  equation  may  be  expressed 


T=  J  © 


=  J 


d© 

dt 


Impulse  is  equal  to  the  time  rate  of  change  of  momentum.  The  ex¬ 
pressions  for  the  rates  of  change  of  momentum  of  the  gears  treated 


separately  may  be  written: 


-i'll, 'A 
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r  j  dt  = 
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(4-30) 


(4-31) 
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Fig.  4- 2a  Schematic  diagram  of  motor  and  load 
inertias  with  interconnecting  gears. 


Fig.  4- 2b  Schematic  diagram  of  mechanical  system 
equivalent  to  Fig,  4-2a. 
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Fig.  4- 2c  Gear  teeth  and  loads  involved 
at  instant  of  impact. 
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where  the  primed  velocities  are  those  following  impact. 

If  the  momentum  is  to  be  conserved  in  the  system  the  impulse 
functions  are  equal  and  opposite 
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r(  ©  - 
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©  ) 
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and 
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from  substitution  of  the  equations 


(4-36) 


rl  =  P2  ,  V  '  P  ®1 

The  assumption  that  all  other  impulse  functicn  s  are  zero  at  the  instant 
of  impact  is  restated  at  this  point. 

A  definition  of  the  coefficient  of  restitution  will  be  made  for  pur¬ 
poses  of  this  work: 
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*  \  ‘i 

p  ©  -  ©  V 
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©  +  ©n 
m  1 


p  ©  -  © 

m  c 


(4-38) 

It  will  be  shown  that  this  definition  of  e  will  satisfy  the  law  of 

conservation  of  energy.  Consider  the  case  of  e=0, 

O  =  —  (p  ©  —  ©!)  ,  p©  =©  (4-39) 

r  m  c  r  c 

in  which  case  for  plastic  impact  (e=0)  the  gears  are  moving  at  the 
same  velocity  following  impact. 

The  expression  for  the  coefficient  of  restitution  is  examined  for 
the  case  of  perfect  elastic  contact,  e  =  l,  for  which  the  law  of  conser¬ 
vation  of  energy  is  satisfied. 
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(4-40) 


Solution  of  the  restitution  equation  for  e=l  yields, 
p  ©_  —  er  =  —  P 

*  JY)  r  m  r> 


?(©  +  ©  )  =  ©  +  © 

r  m  m  c  c 

For  conservation  of  energy 

i  «2  '2  i  . '2  «'2 
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The  momentum  equation  which  is  independent  of  e  is  : 
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_  (0  -©)=(©-©) 
p  J  mm  c  c 

-Lj 

or  rearranging  J 


( ©  -©)  =  (©-©) 
p  J  m  m  c  c 


(4-45) 

(4-46) 


The  factored  energy  equation  is  now  divided  by  the  momentum 
equation  to  yield 


P(©  +©)  =  (©+©) 

»  m  m  c  c 


i 

) 

m  m  e  c 

which  checks  with  the  definition  of  the  coefficient  of  restitution 


(4-47) 


for  e=l . 

The  momentum  equation  implied  from  the  impulse  approach  that 

!  I 

dt  o  hence  ©  =  ©  and  ©  =  ©  •  The  definition  of  the  coefficient 

c  c  mm 

of  restitution  is  independent  of  position.  It  may  be  noted  that  the 
equations  for  energy  and  momentum  depend  only  on  the  inertia  dis¬ 
tribution,  gear  ratio  and  angular  velocities. 
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In  manipulation  of  these  equations  it  is  seen  that  load  inertia 

2 

may  be  transferred  to  the  ,  shaft  by  multiplication  p  ,  the  same 

* 

transfer  may  be  accomplished  with  load  friction.  Since  ©^  ,  and  ©^  , 

are  related  to  ©  and  ©  it  may  be  reasoned  that  the  results  for 
J-j  J _ i 

p=  1  may  be  extrapolated  to  physical  systems  where  p  /  1 

Since  a  total  energy  accounting  may  be  made  for  the  system  with 
the  equations  used,  the  rotational  energy  lost  is  attributed  to  the  heat 
of  deformation  of  the  gear  teeth.  This  transfer  of  energy  could  m 
fact  be  determined  from  the  equations  used.  Thus  the  conservation 
of  rotational  energy  e=  1  is  a  special  case  of  a  broad  interpretation 
of  the  law  of  conservation  of  energy. 
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4 . 3  Computer  Program  Development 

The  physical  equations  of  the  net  systems  were  programmed  for 
solution  using  the  Control  Data  Corporation  1604  high  speed  digital 
computer  utilizing  paper  tape  program  input  and  magnetic  tape  output. 
An  IBM  717  line  printer  was  used  to  extract  data  from  the  magnetic 
tape  The  Control  Data  Corporation  machine  library  and  the  U.  S. 
Naval  Postgraduate  School  computer  subroutine  library  were  used  for 
assembly,  Runge-Kutta- Gill  numerical  integration,  and  decimal 
output 

Several  changes  were  made  in  the  forms  of  the  physical  equations 
of  the  net  system  in  order  to  eliminate  duplicate  computing  operations 
and  to  fit  the  equations  to  a  form  suited  to  the  variable  parameters. 

In  the  table  of  symbols,  computer  mnemonic  (m)  terms  which  are  used 
in  the  assembly  subroutine  have  been  indicated  by  parenthises  and  will 
be  defined  when  encountered. 

The  equation  of  the  system  with  motor  and  load  combined 


(4-48) 


was  put  in  the  form 


where 


(4-49) 
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OMEGANSO  =U)r2 
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(RUNGE)  ,  numerical  integration  method  The  RUNGE  sub¬ 
routine  required  the  definition  of  synonymous  terms  for  the  four 
iterative  cycles  used  to  produce  one  extrapolated  set  of  variables. 
The  increment  of  the  independent  time  variable  chosen  was  0.01  sec. 
for  all^“^0 . 1  and  0 . 004  sec  for  all^* ^  0.1. 


(4-50) 


(4-51) 


The  load  free  equation 

f  L  • 

was  put  in  he  form  ©  +  - —  ©  =  O 

c  J  c 

L 

for  the  computer  0  +  B  @  _  0 

c  —  c 

A  separate  set  of  RUNGE  synonyms  was  used  for  the  load  free  equation 
in  order  not  to  destroy  the  previous  computations  before  they  were 
determined  of  no  further  value  in  computing  and  also  in  the  making  of 

program  decisions 

*  •* 

©  was  designated  V  and  ©^  was  designated  VDOT  for  the  computer 
program  . 

The  equation  for  the  motor  alone 
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m  mm 

was  put  in  the  form 
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(4-54) 


and  for  RUNGE,  the  terms 


©  =  THETAM 

m 

and  ©  =  WDOT  were  used, 
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The  equation  representing  the  law  of  conservation  of  momentum 


and  the  definition  of  coefficient  of  restitution  were  combined  to  the 

2 

.T  .  it 

forms 
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for  the  computer  program  No  provision  for  additional  m  terms  was 

made  to  denote  the  primed  values.  The  additional  terms: 

2 

e  =  RESTITUT  p  =  RHO  and  p  ~  RHOSQ  were  used , 

Two  equations  were  used  to  define  the  boundaries  of  operation  of 

the  load  and  motor  when  they  were  acting  separately, 

©  =  ©  and  ©  =  p  ©  +  A  (4-57) 

c  »  m  c  I  m 

Two  major  decisions  of  the  computing  cycle  were  the  determination 
of  the  point  where  the  motor  and  load  would  float  free,  F LOATEST  , 
and  the  response  of  the  system  to  impact  of  the  load  and  motor  gears 
when  the  boundary  conditions  on  position  were  met,  COM3  TEST  0 

The  first  decision,  the  point  of  float  free,  is  made  on  the  equality 
or  inequality  of  slopes  in  the  phase  plane,  i.e*  ,  at  separation 


N  =  N  Equating  the  slopes: 

L  s  2 

un 

__(OR~0c)  +  B_  =  O 


(4-58) 


The  above  equation  is  solved  after  computing  each  point  in  the 
combined  phase  trajectory.  If  N^*  O  the  system  remains  com¬ 

bined.  If  N  +  N  <  O  the  sytem  separates,  having  the  response 
S  L  "Sh 

motor  alone  and  load  alone  immediately  thereafter. 
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The  second  major  decision,  that  of  whether  after  impact  the 
response  should  be  that  of  the  combined  system  or  the  response  of 
the  load  and  motor  acting  separately,  is  made  on  the  basis  of  the 
resultant  velocities  after  impact  and  upon  which  side  of  the  backlash 
the  motor  and  load  positions  are  found. 

The  computer  program,  its  flow  diagram,  a  description  of  the  pro¬ 
gram  operation,  and  an  explanation  of  the  readout  of  information  with 
examples  are  presented  in  the  Appendix, 
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4 . 4  Methods  and  Scope  of  Examination 


Two  modes  of  operation  of  the  computer  program  were  used 
The  first  mode  considered  was  the  phase  trajectory.  This  mode  was 
used  to  check  the  operation  of  the  computing  cycle.  Time,  velocity 
and  position  of  the  load  were  printed  when  the  system  was  combined 
When  the  system  was  separated,  the  values  of  motor  velocity  and  motor 
position  were  also  printed.  Although  the  time  increment  used  for  com¬ 
putations  was  either  0,01  sec.  or  0.004  sec.  ,  printed  outputs  for  the 
first  mode  were  taken  only  every  0. 1  sec.  Additional  printed  outputs 
were  taken  at  the  time  of  contact  of  the  gears  just  prior  to  and 
immediately  following  the  solution  of  the  momentum  and  restititiL;.'* 
equations.  A  typical  printout  of  the  phase  trajectory  mode  of  operation 
is  shown  in  the  Appendix . 

The  second  mode  of  operation  was  used  for  priv  Ang  only  the  maxi¬ 
mum  computed  positive  load  overshoot  position  for  each  cycle,  the 
associated  load  velocity  and  the  exact  problem  time  of  the  computation. 
A  sufficient  number  of  print- outs  for  each  problem  was  obtained  to 
prove  limit  cycle  existence  and  average  size  Since  this  investigation 
was  concerned  primarily  with  the  steady  state  response  the  second 
mode  of  operation  was  the  one  utilized  to  obtain  the  majority  of  the  data. 
This  read-out  method  markedly  decreased  the  data  reduction  time  for 
the  problem  analysis,  as  opposed  to  the  analysis  of  steady  conditions 
provided  by  a  full  phase  trajectory  print-out  Since  the  time  required 
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for  computer  read-out  was  several  magnitudes  greater  than  the 
computing  time,  valuable  computer  time  was  saved  by  this  mode  of 
operation.  An  example  of  the  printed  output  obtained  when  only  maxi¬ 
mum  load  position  overshoot  was  of  interest  is  given  in  the  Ajoendix. 
Approximately  1400  phase  trajectories  were  solved  using  this  mode  of 
operation.  Each  solution  required  an  average  of  three  minutes  of  com¬ 
puter  time. 

For  the  purpose  of  general  examination  of  the  problem  it  was 

2 


assumed  that  p  =  1,  K=  l,U)n=  1,  p  J^=  lr  Specific  solutions 

were  made  with  p  ^  1,  1  to  determine  scaling  effects.  It  was 

determined  that  velocities  were  scaled  in  direct  proportion  with  uJ^ 

and  time  is  scaled  in  an  inverse  ratio  with  u)  It  was  also  determined 
2  fL  2  JL 

that  p  ^ —  and  p  — —  were  valid  parameters  and  in  agreement  with 
*  T  m 

the  non  dimensional  parameters  used  for  the  investigation. 

The  program  can  also  be  made  to  solve  linear  systems  by  setting 


=  O. 


The  major  variable  parameters  used  in  this  investigation  were: 


^  =  0.1,  0.2,  0,3,  0.4,  0.5,  0.6,  0  8,  1.0 
2f 

P  f  T 

— —  =  0  ,  0.2,  0.4.  0,6,  0.8,  1.0 

ft 


2 

P 


m 


=  0.1. 
0.9, 


0.2,  0,5,  0,8,  0.9 
0.8,  0.5,  0.2,  0.1 


=  0 . 3  radians 


e  =0,0.6,08,10 
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Additional  parameters  of 


p2< 


=  0 . 04,  0,1,  0 . 9;  and  ^=  0  01,0.03,  0,1,0.15  radians  were 


used  in  certain  instances  to  examine  particular  characteristics  of  the 
response . 

The  program  was  limited  somewhat  in  that  the  value  of 


y—  =Oo  and  zero  were  excluded  due  to  generation  of  undefined 
m 

mathematical  quantities. 

The  value  of  backlash  was  made  abnormally  large, ^  =0.3  radians, 
to  allow  an  easier  interpretation  of  the  non-linear  response.  S-ince 
the  influence  of^was  linear,  this  caused  no  inaccuracies.  All 
graphs  for  limit  cycle  size  were  plotted  with  ^  =  0,3  radians.  The 
effect  of  values  for  other  than  0, 3  is  also  shown  in  graphical  form. 
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5.  ANALYSIS  CF  AFSULTS 

The  primary  purpose  of  this  study  was  the  determination  of  the 
conditions  for  existence  of  a  limit  cycle.  Due  to  the  nature  of  the 
computations  an  analysis  of  the  results  of  each  numeric  case  pro¬ 
vided  an  individual  nyesM  or  "no"  answer,  i  e  ,  a  limit  cycle  existed 
or  did  not  exist.  This  information  plots  conveniently  as  individual 
points,  and  has  been  recorded  on  Figs.  5-1  and  5-2  in  condensed, 
semi- quantitative  form  to  give  a  broad  view  of  limit  cycle  conditions. 
On  Figs  5-1  and  5-2  the  small  circles  merely  indicate  points  at 
which  a  computation  was  made,  the  curved  lines  were  averaged  in 
as  accurately  as  possible  with  the  available  data.  In  general  these 
figures  are  sufficiently  accurate  for  an  engineering  check,  since 
parameter  values  which  are  close  to  a  boundary  obviously  lead  to 
a  long- duration  oscillation  even  if  there  is  no  limit  cycle. 

The  data  in  Figs.  5-1  and  5-2  can  be  presented  in  the  form  of  a 
three  dimensional  picture  which  provides  considerable  perspective. 
Such  a  presentation  is  given  on  Fig.  5-3,  which  is  a  three  dimensional 
sketch  of  the  limit  cycle  boundary  for  the  particular  case  of  e=  1 . 0 
The  surface  of  this  contour  is  the  limit  cycle  boundary.  For  com¬ 
binations  of  parameters  specifying  a  point  on  or  below  the  surface 
there  is  no  limit  cycle.  For  combinations  of  parameters  specifying 
a  point  above  the  surface  a  limit  cycle  is  guaranteed,  and  the 
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Limit  Cycle.  Contour  Chdrll 

(not  to  SCRt-E.  SEE  PtCr.  5-&) 

CijclC  —  O  radians 


amplitude  of  the  limit  cycle  is  related  to  the  distance  from  the  surface 
to  the  pertinent  point.  This  is  further  elaborated  on  by  the  sketches 
of  Fig.  5-4. 

For  calculation  purposes  detailed  and  quantitative  curves  are 
desired,  and  the  three  dimensional  figure  is  not  a  convenient  tool. 
Therefore  the  results  of  the  digital  computer  study  are  summarized 
using  a  series  of  families  of  curves.  These  are  given  in  Figs.  5-5 
through  5-28  to  show  the  effects  of  parameter  variation  on  both  the 
existence  and  amplitude  of  the  limit  cycle.  This  is  essentially  the 
data  which  was  summarized  in  Figs.  5-1  and  5-2  (limit  cycle  exist¬ 
ence),  but  is  more  complete  and  more  accurate. 

In  all  of  the  computations  reported  to  this  point  the  backlash  was 
set  a^/\  =0. 3  radians.  It  was  tacitly  assumed  that  all  amplitude  pheno¬ 
mena  would  be  directly  proportional  to  the  amount  of  backlash.  Since 
this  condition  was  easily  checked  by  additional  computations  it  was 

decided  to  verify  the  assumption.  The  results  are  shown  on  Fig.  5-29 

# 

and  5  30,  from  which  it  is  evident  that  the  amplitude  of  the  limit  cycle 
is  indeed  directly  proportional  to  the  amount  of  backlash. 

Included  for  purposes  of  general  information  are  two  phase  tra¬ 
jectory  plots  shown  on  Fig.  5-31  and  5-32.  These  show  the  type  of 
information  available  from  the  digital  computer  data. 
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Limit  Cycle  Error  (Radians) 


Fig.  5-5 
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Limit  Cycle  Error  (Radians) 


Fig.  5-6 
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Fig.  5-7 
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Limit  Cycle  Error  (Radians) 
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Fig.  5’ 9 
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Limit  Cycle  Error  (Radians) 


Fig.  5-lZ 
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Fig.  5-13 
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Limit  Cycle  Error  (Radians) 


Fig.  5-15 
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6.  APPLICATION  OF  RESULTS 

In  this  section,  the  results  are  applied,  to  hypothetical  physical 
system s .  Methods  are  indicated  by  which  physical  systems  may  be 
analyzed  or  designed. 

It  is  desirable  to  find  an  equation  relating  certain  system  variables, 
therefore  the  equations’ 


UJn  - 


and 


p  K  p  K 

2 

J  +  p  J  J 
m  r  L  m 
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1  +  p  J. 


£  +  of 

m  r 


?  £ 
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m 

i 


+ 1 


J  +  P  JT 


m 


^  +  P2j  L 

m  L, 


will  be  mampalated  to  the  form 

f 


m 


2 

P  f 


m 


P2f 


+  1 


'/Jm  s/‘ 


+  P  J 


r  p 


m 


2  _  - 


1  +  p  J 


or 


\/p  vK 


m 


p2f  i 


m 


+  1 


J  K  (l  +  p2  J  \ 
m  ^  r  L  / 


m 


m 


1/2 


(6-1) 


(6-2) 


(6-3) 


(6-4) 


In  view  of  the  form  of  Equation  6-4  for  and  recalling  that  the 

r  2  2 

variable  parameters  used  on  the  figures  were  p  J  and  p  f  = 

L  L 


m 
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it  is  seen  that  the  expression  for 


2 

P 


£ 


L 


f 

m 


+  P 


f 

m 


+  1 


in  fact  includes  the  other  two  variable  parameters. 

The  value  of  overshoot  and  hence  limit  cycle  are  not  affected 
by  Ul0;  however,  the  time  and  velocity  are  scaled  for  UJn /  1.0.  Thus 
the  transient  response  characteristics  may  be  related  to  the  same 
parameters  used  in  these  applications. 

If  a  motor  and  load  are  selected  for  a  system  and  the  limitation 
is  imposed  that  system  operation  shall  not  result  in  a  limit  cycle. 
Equation  6-4  and  the  limit  cycle  existence  charts  may  be  used  to 
determine  the  unspecified  system  parameters  ,  K  and  p.  The 
coefficient  of  restitution  for  the  proposed  gear  train  material  must 
also  be  determined. 

To  proceed,  select  a  value  of  K  and  p,  then  determine  the  value 

of  £  from  Equation  6-4.  Enter  the  limit  cycle  existence  chart  for 

.  2 

the  value  of  e,  system  and  the  parameters  p  and 

T~ 


f 

m 

Tt, 


+  1 


Determine  if  the  selected  parameters  describe  a  point  of  no  limit 
cycle.  Since  limit  cycle  existence  is  independent  of  backlash  O 

excluded),  backlash  is  not  a  variable  parameter.  Successive  trial 
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values  of  the  unspecified  parameters  may  be  required  to  meet  the 
limitation  of  no  limit  cycle.  If  it  is  also  required  that  certain  transient 
response  characteristics  be  met,  these  conditions  cc  ild  be  examined 


by  a  similar  method. 

In  the  event  that  no  acceptable  parameters  are  found  to  satisfy 
the  system  requirements  for  no  limit  cycle,  the  figures  for  limit  cycle 
size  may  be  examined  to  obtain  a  minimum  size  limit  cycle.  The 
method  of  solution  is  similar  to  that  previously  described. 

Example:  The  material  intended  for  use  in  the  gears  has  an 
e=  0.6.  The  motor  and  load  have  the  parameters  f  =  0 . 64,  J  =0.25, 


m 


m 


f  =  0.64,  J  =1.0  and  it  desired  that  the  system  operate  with  no 
L  L 

limit  cycle.  If  the  values  are  selected  K-4.0,  p=  0.5,  the  values 


p  J  T  =1.0  and  f  =4.0, 
1  >  m 


f 


=  0.2  result. 


m 


P2  £1 


m 


P 


+  1 


3/2 


Solve  Equation  6-4  for  =  .  5  (  .  64)  ( 5+  1) 

[25  (4)  (l+lj1/2 

=  0.4.  Enter  the  chart  for  limit  cycle  existence,  Fig.  5-1 
e=0. 6,^*  =  0.4, 


,  thus 


f 


=0.2,  p  J  =  1.0 


m 


P2fi 


+  1 


m 


It  can  be  seen  that  limit  cycle  will  exist  for  these  parameters. 
Since  K  is  most  easily  varied  in  equation  6-4,  a  new  K  is  selected 
at  1 .0  in  an  effort  to  find  parameters  for  no  limit  cycle.  It  is  seen 
that 


C  =  C  A  =  0.4  (2) 

3  K  =  1 .0  3  K=4.0  N 
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with  all  other  parameters  constant,  thus£*  =0.8  Enter  Fig.  5-1 
with  e=  0 . 6 


£=o.sPV=i.o  1 


f 


=  0.2 


m 


m 


P2fi 


+  1 


These  parameters  describe  system  operation  with  no  limit  cycle  and 
are  well  within  the  area  limits.  The  transient  response  of  this  system 
may  be  undesirable  and  various  combinations  of  p  and  K  would  have 
to  be  tested  to  satisfy  the  additional  requirements  of  the  problem.  It 
will  be  restated  that  the  existence  areas  of  Figs.  5-1  and  5-2  are  not 
exact. 

The  figures  of  limit  cycle  size  may  be  used  in  the  same  manner  as 
the  existence  charts.  The  same  parameters  of  the  first  case  are 


assumed  :  k=4.0,  p=0.5,  pj  =  1.0, 

L 


1 


=  0.2 


m 


m 

PZ{ 


+  1 


r=0.4  -  -L 

e=0.6,A  =  0.1  radians.  Enter  Fig.  5- 12  with  the  parameters  given. 
The  ordinate,  Limit  Cycle  Error  (0.3)  =0.13  radians,  is  obtained. 

A 

The  maximum  positive  error  in  the  limit  cycle  =  0.13  (0.1)=  0.0434 

0.3 

radians.  Enter  the  same  figure  for  e=0.6,^*=0.8 


?\  -  '-0 


f 


=  0.2,  A  =0.5  radians;  this  is  the  second 


m 


m 


P2fi 


+  1 


case 


of  examinati  on  f  or  existence.  No  ^  -  0.8  curve  exists;  hence  no 
limit  cycle  exists  for  this  choice  of  parameters  for  any  nominal  value 
of  backlash. 
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7 


•  >  A  C;TIC  IMPACT 


Tr>  A  CHATw  ■  rTEF.KTTJCS,  ~ 

The  effects  of  backlash  on  the  transient  response  to  a  step  displace¬ 
ment  input  is  also  of  interest,  so  a  relatively  exhaustive  study  was  made 
of  the  plastic  impact  case.  The  char  cteristics  of  interest  are  the  peak 
overshoot  Mpt,  the  settling  time  and  the  peak  time.  Of  these,  the  peak 
overshoot,  Mpt,  produced  the  most  satisfactory  pattern  with  the  result 
that  d<  ta  charts  and  equations  have  been  prepared  to  assist  in  analysis 
and  design  Settling  time  and  peak  time  were  much  less  satisfacv  .  y  to 
work  with  The  results  obtained  to  date  are  presented  here,  but  more 
work  remains  to  be  done  . 

The  bulk  of  the  computed  data  assumes  that  the  backlash  is  0. 3 
radians  To  use  this  data  with  other  values  of  backlash  it  is  first  necessary 
to  eetablish  the  effect  of  varying  the  amount  of  backlash.  This  is  done  on 
Figs.  7-1  *2,3.  These  figures  merely  verify  the  intuitively  obvious 
fact  that  the  amount  of  overshoot  must  vary  linearly  with  the  backlash, 
but  that  the  proportionality  constant  depends  on  the  various  parameter 
values 

7.1  Maximum  Overshoot;  Plastic  Impact  Case 

Figs.  7-4, 5, 6,7, 8,  and  9  present  the  results  of  computer  calcu¬ 
lations  in  the  form  of  families  of  curves.  Each  family  is  a  plot  of 
2  2  f  * 

Mpt  vs  p  J  /  J  for  a  specific  p  /FT  ratio  with  $  as  the  family 

LI  m 


K  &  A  Figs.  5,6,7 

Figs.  10, 11 , 12, 13, 14, 15  in  K  &  A 
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parameter.  All  are  for  the  case  of  plastic  impact  These  curves  can 
be  used  to  predict  the  peak  overshoot  of  any  arbitary  second  order  system 
with  backlash  and  plastic  impact  by  reading  off  the  pertinent  values  and 
inserting  in  the  equation 


Mpt 


A=Ai  Mpt 

A=0  +  0.3 

Mpt 

A=  .3  ~  MPt 

srssi 

£  ~  S’  i 

A=D 

S=Si 


(7-1) 


As  was  pointed  out  earlier,  the  solutions  were  made  and  the  results 

were  plotted  using  a  backlash  angle  of  0  3  radian.  Thus  the  values  of 

amplitude  of  maximum  overshoot  on  the  curves  are  significant  only  in  that 

they  represent  the  slope  of  the  overshoot  versus  backlash  straight  lines. 

This  slope  must  be  calculated  in  order  to  scale  the  value  of  overshoot  read 

from  the  curve  for  the  magnitude  of  backlash  actually  expected  Therefore, 

in  order  to  refine  the  prediction  scheme,  other  design  curves,  figures  7-10 

through  7*15,  were  calculated  using  the  curves  of  figures  7-4  through  7-9 

dM 

just  discussed.  Along  the  curves  the  slope  pt  was  calculated  by, 

d  A 


dM 

pt 


dA 


M 

pt 


4=0 


i 


.  3 


(7-2) 


These  values  of  dM  were  then  plotted  versus  inertia  ratios  at 

Pt 

constant  friction  ratios  before.  These  curves  give  a  more  straight¬ 
forward  prediction  method  with  the  following  equation' 


K  &  A  ,  16-21 
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99 


M 


Pt 


=  M 


A=AJ 

r*r 


pt 


A=o 

r=?, 


+  dM  /  d  A| 
pt 


k-ct 


x  A 


(7-3) 


In  order  to  demonstarte  the  procedure  for  using  the  design  charts 
shown  in  figures  7  4  through  7-15,  &n  example  problem  is  solved  obtaining  the 
maximum  overshoot,  M  ,  for  a  second  order  servo  with  backlash.  The 

Pt 

differential  equation  for  the  combined  motion  of  an  armature  controlled 
motor  is: 


©  + 
c 


f  +  f 
t  1 


j  +p  j 

m  '  l 


©  + 
c 


PK  a 


J  J  T 

m  Lj 


pK  a 

©  =  -L-±  e  (7-4) 

L  _  c  T  .  r 


J  +P  J 
m  r  L 


Given  the  following  specifications  of  a  D,  C  „  shunt  motor  coupled 
through  a  gear  box  to  a  load: 

Horsepower  =  .008  ? 

Rated  RPM  =  4000 


R 

282  ohms 

K 

=  .233  volts/ radian/ second 

g 

J 

-6  2 
=  8.2  X  10  slug  feet 

m 

f 

=  3.  34  X  10‘6  lb. -feet/ radian/ 

m 

K 

=  10  volts/  radian 

e 

P 

=  i 

A 

=  05  radians  1 

j 

-6  r  2 
=  16.4X10  slug-feet 

L 

f 

=  10.4X10  lb  . -feet/ radian/ 

L 

K 

=  .  1464  lb. -feet/ampere 
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The  required  values 


K  = 
a 


K  K 
t  e 

R 


shown  in  Equation  7-4  are  obtained 

■  1464  x  10  _  .0052  LB-FT 

282 


as  follows- 


14  5 


RAD 

SEC 


13.74  x  10 


LB-FT 

RAD/SEC 


ft=  fT  +  fi  =  <13-74  +  121)  x  10 


135  x  10"6 


LB-FT 

RAD/SEC 


K  K 

f  _  t  g  _  _  1464  x  .233  _  121  x  10 

1  R  282 


LB-FT 

RAD/SEC 


2 

P 


JT  /  J 

1^  m 


(lp  1 6 . 4  x  10 


8.2  x  10 


Tr 


2.0 


10.4  x  10  ~6  _  .077 

135  x  10  '6 


Two  methods  are  shown  for  obtaining  the  maximum  overshoot  with  a 
backlash  angle  of  .05  radians  in  the  gear  trains. 

The  first  method  is  accomplished  using  the  curves  of  figures  7-4 
through  7-9  The  procedure  is  simple,  but  requires  more  steps  than  does 
the  second.  Also  one  must  remember  that  these  curves  are  drawn  for  a 
backlash  of  0  3  radian  and  must  be  scaled  for  the  expected  amount  of  back¬ 
lash.  The  second  method  makes  use  of  the  design  charts  of  figures  7-10 
through  7  -  15  . 
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Method  1 


Z  Z 

Entering  figure  7  4  (  n  f  /  F_  =  0  1) ,  at  o  JT/J  =  2.0,  and  a 

1  L  /  T  (  L  m 

damping  ratio,  =  0.19,  the  amount  of  overshoot  for  a  backlash  angle,  £  = 


0 . 3  rad  .  is  M  =  1 . 658  .  Then: 

pt 


M 


pt 


=  Mpt 

+  .05/. 3 

M 

pt 

M 

pt 

£  =  •05 

II 

o 

£=.3 

r=!9 

r=  19 

o 

II 

A=o 

r=  ■  19 


From  second  order  linear  characteristics 
M 


hence , 


pt 


pt 


Method  2 


M 


pt 


=  1.545 

=1.545  + 

.05/ . 3 

1.658  -  1 . 545 

=  1 . 545  + 

.0189 

=  1.564 

=  M 

+  dM  /dA 

.£=•05 

pt 

A  =0 

r=-i9 

f  =  .19 

X  A 


2  2 

from  figure  7  10,(p  £  J  FT  =  0 . 1) ,  at  p  J  J Jm  =  2  •  0  -  and  a  dampinj 

ratio  J"*  =0.19,  dM  /  d£  =  0.357 

pt 


Also: 


M 


pt 


=  1.545 


A=0 
f  =.19 


hence , 


M  =  1 . 545  +  .  357  X  .  05 

pt 

=  1.545+  .018 
=  1 . 563 
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The  procedure  for  using  the  charts  rem?;ns  unchanged  for  servo 
systems  which  utilize  other  types  of  motors.  The  only  difference  which 
arises  is  in  the  calculation  of  the  chart  parameters  from  the  known 
system  constants.  For  a  field  control  motor  if  the  field  time  constant 
is  negligible,  the  equation  is  similar  to  the  armature  control  motor  equation. 
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(  7-5) 


In  like  manner  the  equation  for  the  system  with  a  two  phase  motor  is 


©  + 
c 


F  K 
t  n 


J  +  p  J 


m 


r  JL 


©  + 
c 


p 


K 


m 


©  = 


PK 


m 


J  +p  J 


2  _  c  -  .  2 


(7-6) 
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where 


K  =  - 
n 

dT 

dn 

(from 

characteristic 

curves  of  two  phase  motor) 

K  =  ■ 

m 

(from 

characteristic 

curves  of  two  phase  motor) 

<3e  1 

The  calculation  of  the  chart  parameters  is  a  simple  matter  once  the  defin¬ 


ing  equations  are  known. 
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7  2  Maximum  Overshoot,  Elastic  Impact  Case. 

Investigation  of  a  limited  number  of  cases  reveal  that  the  amount  of 
overshoot  with  plastic  impact  is  always  greater  than  that  with  elastic 
impact.  The  difference  between  this  maximum  overshoot  decreases  as 


me rease s  . 


However,  under  no  circumstances  was  the  effect 


onmaximum  overshoot  greater  than  .01  radians.  The  effect  on  the  time 
of  maximum  overshoot  appears  negligible.  In  a  few  cases  the  time  was 
different  for  the  plastic  and  the  elastic  cases  but  the  difference  was  very 
small  and  no  pattern  of  difference  could  be  ascertained. 

These  results  point  out  the  validity  of  the  original  assumption  made 
for  this  study.  In  a  practical  system  there  will  exist,  of  course,  some 
amount  of  gear  bounce  when  the  gears  contact  each  other.  The  effect  on 
maximum  overshoot,  however,  is  so  small  that  the  curves  obtained  m 
this  study  remain  reasonably  accurate 
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7.3  Settling  Time 


In  determining  the  effect  of  backlash  on  the  settling  time  of  a  second 
order  servo  system,  it  was  necessary  to  first  find  a  definition  for  settling 
time  which  would  be  suitable  for  all  of  the  cases  which  were  investigated. 
Normally,  settling  time  is  defined  as  the  time  for  the  output  to  reach  98% 
of  it  final  value.  This  definition  has  real  meaning,  however,  only  for 
linear  stable  systems  in  which  there  are  no  sustained  oscillations.  For 
systems  which  do  exhibit  small  sustained  oscillations,  the  settling  time  is 
sometimes  defined  as  the  time  required  for  the  amplitude  of  oscillations  to 
decrease  to  within  2%  of  the  final  value.  For  very  small  amplitudes  of 
oscillations  this  definition  gives  an  unrealistic  value  for  settling  time. 

For  the  purpose  of  this  study  settling  time  is  considered  as  the  time  re¬ 
quired  for  the  system  output  to  settle  to  98%  of  the  difference  between  the 
magnitude  of  the  input  and  the  amplitude  of  the  sustained  oscillations. 

For  all  cases  for  which  the  damping  coefficient  was  less  than  0.6 
there  were  no  conclusions  as  to  effects  of  backlash  on  the  system  settling 
time  since  the  magnitudes  of  the  sustained  oscillations  were  always  quite 
large  for  the  amount  of  backlash  used  in  the  solutions  and  the  settling  time 
was  difficult  to  obtain  accurately.  For  those  cases  for  which  the  damping 
coefficient  was  greater  than  0.6,  the  sustained  oscillations  were  small 
if  they  existed  and  an  analysis  revealed  the  following  general  conclusions: 
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1*  As  the  size  of  the  backlash  angle  increases,  the  settling  time 
will  increase  provided  the  other  parameters  are  kept  constant.  For  the 
cases  investigated,  this  increase  in  settling  time  was  nearly  proportional 
to  the  backlash  angle. 

2 

2  As  the  ratio  n  JT  /J  increases,  with  other  parameters  constant, 
*  L  m  r 

the  settling  time  increases. 

2 

3.  At  the  ratio  *  f  /  F  is  increased,  with  other  parameters 
v  L  t 

constant,  the  settling  time  decreases. 

The  data  included  in  Tables  I  -  III  show  these  general  conclusions 
and  also  some  quantitative  effects  of  backlash  on  settling  time. 
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TABLE  I 


SETTLING  TIME  DATA 

Hr1-0 


f 

“M 

p2fL 

TO  ' 

Radians 

Seconds 

S 

.6 

1/9 

.2 

.3 

28 

.6 

1/9 

.4 

.3 

21 

.6 

1/9 

.6 

.3 

16 

.6 

1/9 

.8 

.3 

13 

.6 

1/9 

1.0 

.3  - 

7 

.6 

1A 

.2 

.3 

28 

06 

i/4 

.4 

.3 

24 

.6 

iA 

.6 

.3 

16 

.6 

1/4 

.8 

®3 

14 

.6 

1/4 

1.0 

.3 

9 

.6 

1 

.2 

.3 

29 

.6 

1 

.4 

.3- 

28 

.6 

1 

.6 

.3 

23 

.6 

1 

.8 

.3 

17 

.6 

1 

1.0 

.3 

10 

.6 

4 

.2 

.3 

29 

.6 

4 

.4 

.3 

28 

.6 

4 

.6 

.3 

26 

.6 

4 

.8 

.•3 

23 

.6 

4 

1.0 

.3 

10 

.6 

9 

.2 

.3 

29 

.6 

9 

.4 

.3 

29 

.6 

9 

.6 

.3 

26 

.6 

9 

.8 

.3 

23 

.6 

9 

1.0 

.3 

10 

TABLE  II 


SETTLING  TIME  DATA 


.  - 1.*  0 


f 

/jM 

FT 

Radians 

A 

Seconds 

t 

& 

.6 

1/9 

.8 

.03 

5 

.6 

1/9 

.8 

.10 

7 

.6 

1/9 

.8 

.15 

8 

.6 

1/9 

.8 

.30 

13 

.6 

1 

.8 

.03 

11 

.6 

-i 

JL  » 

oQ 

.10 

14 

.6  ' 

n 

± 

.8 

.15 

16 

.6 

1 

.8' 

.30 

17 

.6 

4 

.8 

.  .03 

12 

.6 

4 

.8 

.10 

14 

.6 

4 

.8 

.15 

17 

«6 

4 

.8 

.30 

23 

.6 

9 

.8 

..03 

13 

.6 

9 

.8 

.10 

16 

.6 

9 

.8 

.15 

18 

."6 

9 

.8  . 

.30 

23 
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TABLE  III 


SETTLING  TIME  DATA 
LJr»  =1*0 


y 

p  JL/ 

/jM 

FT 

Radians 

Seconds 

t 

B 

.8 

1/9 

.2 

.3 

8 

.8 

1/9 

.4 

.3 

.7 

.8 

1/9 

.6 

.3 

6 

.8 

1/9 

.8 

.3 

5 

.8 

1/9 

1.0 

.3  * 

5  . 

.8 

'  1/4 

.2 

.3 

20 

.8 

1/4 

•4 

.3 

9 

.8 

1/4 

.6 

.3 

7 

.8 

1/4 

.8 

.3 

6 

.8 

1/4 

1.0 

.3 

5 

.8 

1 

.2 

.3 

43 

.8 

1 

.4 

.3 

31 

.8 

1 

.6 

.3 

10 

.8 

1 

.8 

.3 

8 

.8 

1 

1.0 

.3 

6 

.8 

4 

.2 

.3 

48 

.8 

4 

.4 

.3 

36 

.8 

4 

.6 

.3.. 

12 

.8 

4 

.8 

.3 

9 

.8 

4 

1.0 

.3 

6 

.8 

9 

.2 

.3 

32 

.8 

9 

.4 

.3 

43 

.8  . 

9 

.6 

.3 

13 

.8 

9 

.8 

.3 

9 

.8 

9 

1.0 

.3 

6 

y-ig  v 
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7.4  Time  of  Maximum  Overshoot 


For  each  of  the  cases  investigated  the  time  of  maximum  overshoot 
(peak  time)  was  recorded  from  the  printout  of  the  solution.  It  was 
found  that  the  presence  of  backlash  causes  a  slight  increase  in  the  peak 
time  over  the  linear  system  for  those  cases  where  .  5  ,  and  has 

negligible  effect  in  those  where  J  <^.5.  When  p^f^  ls  increased  from 


ze 


ro  to  0  5,  other  parameters  constant  and^** ^  5  ^  peak  time  increases 


slightly  and  remains  essentially  constant  for  p  f  ^>5.  As 


h 

J  /J  is  increased,  other  parameter s  constant  and  (*  .  5  ,  peak  tim< 

T  m  ^ 


PJL 

is  decreased. 
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7.5  Other  System  Characteristics. 


A  careful  study  of  the  results  as  plotted  in  figures  7-4  through  7-9 
reveals  some  interesting  phenomena.  At  all  friction  ratios  it  is  observed 
that  as  the  ratio  of  load  inertia  to  motor  inertia  is  made  smaller  the  ampli¬ 
tude  of  the  overshoot  is  less.  There  is  a  point  where  the  overshoot  is  identi¬ 
cal  to  that  of  a  linear  second  order  system,  in  other  words,  where  backlash 
has  no  effect  on  the  overshoot  Note  that  this  point  occurs  at  different 
inertia  ratios  for  each  value  of  damping  ratio.  Also  this  point  varies  de¬ 
pending  on  the  value  of  friction  ratio  If  the  inertia  ratio  is  decreased  beyond 
this  point  then  the  amplitude  of  overshoot  remains  constant  at  the  value  of 
the  linear  case  Figure  7-16  is  a  plot  of  the  boundaries  of  the  combinations 
of  variables  which  will  cause  an  effect  on  the  amplitude  of  overshoot  if  the 
system  has  any  appreciable  amount  of  backlash  The  coordinates  are  fric¬ 
tion  ratio  and  inertia  ratio  with  the  boundaries  labeled  with  values  of  damping 
ratio.  Areas  to  the  right  of  a  boundary  are  combinations  of  the  variables 
which  affect  overshoot.  These  curves  were  plotted  using  information  deducted 
from  figures  7-4  through  7-9.  Three  points  were  selected  as  test  points  to 
prove  the  validity  of  the  boundaries  as  plotted.  These  points  are  shown  on 
the  figures  as  points  nAM  ,  ,!Bn,  and  "CM.  The  system  parameters  used  at 
the  test  points  were: 

Point  A  Point  B  Point  C 


F 


T 


=  0.8 


pJjJ  =0.25 

p  L  m 


1  14 


-t- 

~t- 

jtif 

in 

p~D.T 

I 

~T 

n3 

nil 

1 

1 

The  system  was  solved  on  the  digital  computer  to  obtain  the  maximum 
overshoot  for  these  three  sets  of  parameters.  According  to  the  curves 
point  A  falls  to  the  left  of  the  boundary  of  the  r  =  .8  area.  Therefore, 
the  overshoe  t  should  be  the  same  as  that  of  the  linear  case.  The  com¬ 
puter  solution  gave  an  overshoot  of  1.019,  which  is,  in  fact,  identical 
to  that  of  the  linear  case.  Similarly,  point  B  falls  to  the  right  of  the 

=0.6  boundary  and  in  the  area  where  the  amplitude  of  overshoot  should 
not  be  the  same  as  for  the  linear  case.  The  computer  solution  gave  a  maxi¬ 
mum  overshoot  of  1.110  which  is  slightly  greater  than  the  linear  case. 

Point  C  falls  to  the  right  of  the  j  =0.3  boundary  and  in  the  area  where 
the  amplitude  of  overshoot  should  be  different  from  that  of  the  linear 
system.  The  computer  solution  gave  a  maximum  overshoot  of  1.41  while 
that  of  the  linear  system  is  1.372.  From  these  boundary  curves  one  can¬ 
not  determine  the  exact  effect  of  backlash  on  the  maximum  overshoot  of  a 
system,  but  can  determine  in  which  systems  backlash  will  cause  an  appreci¬ 
able  effect. 

Another  phenomena  which  is  revealed  in  the  results  as  plotted  in 
figures  7-4  through  7-9  is  that  for  all  cases  the  greatest  maximum  over¬ 
shoot  occurs  in  the  range  0.1^ p^J  ^  10.0  for  the  range  of 

2 

a  J  /J  tested.  However,  it  appears  that  since  the  curves  are  trending 
f  L  m 

2 

upwards  at  the  high  end  of  the  p  J^/J^  axis  that  at  very  large  inertia 
ratios  the  overshoot  would  be  greater.  It  is  felt  that  these  large  inertia 


ratios  are  outside  the  practical  range.  An  attempt  was  made  to  explain 
these  phenomena  but  without  success.  It  is  often  very  difficult,  if  not 
impossible,  to  find  a  physical  explanation  for  certain  behavior  caused  by 
non-linearities  They  can  merely  be  pointed  out  from  analyses  of  experi  - 
mental  information  and  accounted  for  accordingly 

It  is  interesting  to  note  that  when  the  ratio  of  load  friction  to  motor 
friction  is  below  0.4,  as  m  figures  7-4  and  7-5  and  11,  there  is  a  possi¬ 
bility  of  an  overshoot  for  the  critically  damped  case,  i.e.  =1  0, 
depending  on  the  inertia  ratio.  At  first  glance  one  would  reason  that  if 
the  system  is  critically  damped  there  should  be  no  overshoot.  However, 
it  can  be  seen  that  if  the  friction  ratio  is  small  and  the  inertia  ratio  is 

large  then  the  ratio  of  F  /J  ,  which  is  the  slope  of  the  drifting  load 

I-*  L 

trajectory,  is  made  small.  Therefore,  when  the  load  separates  from  the 
motor,  its  velocity  remains  fairly  constant  at  a  relatively  high  value  and 
its  momentum  carries  it  past  the  correspondence  point. 
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8.  CONCLUSIONS 


Since  this  report  is  to  be  the  first  of  a  series,  it  is  apparent  that 
the  results  are  incomplete  and  the  conclusions  must  therefore  be  brief: 

1  .  It  has  been  established  by  calculation  the  t  a  second  order  feedback 
control  system  with  backlash  may  or  may  not  exhibit  a  limit  cycle. 

2.  The  existence  of  a  limit  cycle  as  a  function  of  various  parameters 
has  been  established,  and  a  series  of  charts  suitable  for  analysis  and 
design  have  been  prepared. 

3.  The  amplitude  of  the  limit  cycle  depends  on  the  amount  of  backlash, 
and  charts  illustrating  this  dependence  have  been  prepared. 

4.  Peak  overshoot  depends  on  the  system  parameters  and  the  amount  of 
backlash.  Charts  and  equLtions  suitable  for  analysis  and  design  have 
been  prepared 

5.  The  effect  of  elastic  impact  at  the  gear  teeth  is  significant,  but  much 
less  so  than  was  anticipated.  Most  of  the  design  charts  include  this 


effect . 
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APPENDIX  A 

COMP  U  TA  T /O/V  FLOW  O/ A  OP  AM 

for 

CDC  /GO 4  D/G/TAL  COMPUTER 


120 


i 


frvm  MOTOR.  33 


I 

,  121 

•;  '  •  I 


from  S~9,  CO 


front  £7t  43 


UM»TCYC 


Yea 


to 

PXepiAttZ 

69 


to  67 


\  to  COMBINED.  A 4 


from  S£,6 


OH 

fr&OAROU* 4 

■ 

STOP 


yes 


yes 

r° 

63 

exit'  1 


aTOPAVe 


95 


96 

I 


97 


J 

r 

H 

r 

r 

JL 

90 

91 

V 

J 

L1MPRIHT 

9£ 
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r 

93 

64 


return  to  /orOfrom 


YE5  j  j^O 


96 
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APPENDIX  A  (continued) 


EXPLANATION  OF  FLOW  DIAGRAM  BLOCK  SYMBOLS 

,1.  Clear  fault  stop. 

2.  Print  last  overshoot  stored  in  LIMITBUF  (symbol  888). 

Print  out  position  in  phase  trajectory  of  .fault  stop 
(symbol  666). 

5.  End  of  file. 

!<-.  Actuate  fault  stop. 

5.  Load  zeta  from  Index  1  (b^). 

6.  Print  last  overshoot ' stored  on  LIMITBUF  (symbol  888) . 

Print  position  in  phase  trajectory  where  it  was  manually 
stopped  (symbol  1961 ) . 

7.  Stop  if  Zeta  =  0. 

8.  2jr  u»  .  a 

9.  Is  zeta  =  0.1? 

10.  Is  zeta  less  than  0.1? 

11.  Set  up  computation  time  of  0.01  seconds  for  problem  and 
set  up  for  print  at  every  0.1  seconds  when  zeta  is  greater 
than  0.1. 

12.  Set  up  computation  time  of  0.004  seconds  for  problem  and  set 
up  for  print  at  every  0.1  seconds  when  zeta  is  equal  to  or 
less  than  0.1. 

13.  Set  up  limit  cycle  print-out  routine  for  20  overshoots  and 
average  the  last  8. 

14.  Set  up  number  of  limit  cycle  print-outs  and  averaging  routine 
depending  upon  the  size  of  zeta: 
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f 

Overshoots 

Start 

Averaging 

Divide 
by  * 

Comp. 

Time 

^  0.8 

6 

k 

2.0 

0.01 

0.6 

7 

5. 

2.0 

0.01 

0.5 

8 

5 

5.0 

0.01 

0.4 

10 

7 

3.0 

0.01 

0.3 

14 

10 

4.0 

0.01 

0.2 

20 

12 

8.0 

0.01 

0.1 

17 

12 

5*0 

0.004 

*The  last  row  printed  has  this  number  printed  on  the  far  left 
if  the  problem  has  a  normal  solution.  This  number  indicates 
how  many  overshoots  were  averaged.  The  average  value  has  1.0 
subtracted  from  it  and  is  printed  in  the  far  right  of  the  same 
last  row. 

15.  Page  reject. 

16.  Clear  buffers  and  set  up  for  new  problem. 

17.  Restitution,  e,  from  Index  3,(b,). 

3 

Backlash,  ^  ,  from  Index  4  (b^). 

J  m  from  Index  5  (b_ ) . 

5 

Jr  *  Jm  ,  K  Jlir* 

/’*  =JLi  “  ;  ~ 

from  Index  6  (bg). 

(fyfr)( Zfun)(j7)  =  fi.  i  =  & 

-^tt  -/°2fL  -f"  >* 

18.  .  Have  all  values  been  cycled  through? 

19.  Set  STOPNOW  after  this  run  to  stop  absolutely. 
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20 


=  0? 


.  Backlash,  ^  , 


21.  Print  consecutive  number  of  this  parameter  run  and  all  the 


parameters  used: 


A 

V* 

e 

A 

J  m 

Jl 

Jm  +/*£Jl 

fm 

fu 

/-2 

K 

VJn?c 

£>» 

f 

22.  Set  program  so  always  acts  as  combined  system. 

25.  Print  consecutive  number  of  this  parameter  run  and  selected 
parameters  for  linear  syston.: 


2je,6)w=A  Jl  K 


2k.  Solve  the  second  order  equation  by  Runge-Kutta-Gill  numerical 

integration  every  0.01  seconds  problem  time  when  zeta  is  greater 

than  0.1,  otherwise  every  0.004  seconds.  Initial  conditions 

•  •  • 

are  &L  =  1*0,  61-  G,  =  0  and  (0  =  1.0  or  the  pre- 

R  C  C  C 

viously  computed  point.  Next  point  is  computed  by  four 
iterations  of: 

0  c  ~  ^ h  (*©r  ~  0c)  ”  &  j ^  CO  h  0C 

25. ’  Set  EXIT  from  print  routine. 

* 

26.  Load  time,  an^-  Prin^  "buffer. 

27.  Go  to  LIMITCYC  to  print  only  the  overshoot.  Otherwise  print 
the  phase  trajectory  points  for  each  0.1  seconds  up  to  a  maxi¬ 
mum  of  12  cycles. 

m 

28.  Store  0  for  and  ©  m 
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29. 

30. 

31. ' 

32. 

33. 

3^. 


35. 

36. 


37. 


38. 


39. 

40. 

41. 

42. 

43. 

44. 
45- 
4 6. 


Is  problem  time  0.1  sec.? 

Is  problem  time  greater  than  0.1  sec.? 

Print  the  number  of  point  computed,  time,  ©^  and  €>c  . 
EXIT  (automatically  set  for  desired  Jump-out.) 

Store  time,  ©c  ,  and  for  motor  and  load  initial 

conditions. 

Is  slope  =?  Ng? 


*>„+  ik.  . 
,0C  J  L 

Is  (5>  negative? 

£ 

Store  ©c-  A  in  0  m 


-  positive  value? 


Solve 


©c= 


for  e«  and  ©£ 


by  Runge-Kutta-Gill. 


for  and  ©^by  Runge-Kutta-Gill 

0L  ©m  <  o  ? 

©  L  —  A  =  negative? 

Set  EXIT  from  print  routine. 

Set  EXIT  from  print  routine. 

* 

Load  time,  ©e  and  0c  • 

Load  0^  and  (0m  « 

Print  number  of  point  computed  and  stored  values. 
EXIT  (automatically  set  for  desired  Jump-out). 


\ 
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k'J.  Solve 

6  / 

©e  ' 


*/ 

© 


-is-j-  r^im(i  +  e)  +  0c(/«2-A 
Jw+/>*Jl  Y.  Jn> 

Q  -  e(/g<3r*  +  ®c) 


“e) 


rv» 


OO  / 


"Pm 

J  *■« 


O/ 

©m  + 


-y—  (©«.  "  ©<) 

JlVI 


48. 

49. 

50. 

51. 

52. 
55. 

54. 

55. 

56. 

57. 

58. 

59. 

60. 

61. 

62. 
65. 

64. 

65. 

66. 


Set  EXIT  from  print  routine. 

0 

Is  ©c  positive? ' 

Set  up  to  count  one  phase  trajectory  when  ©«  is  negative. 

16  9c  >  1-0? 

Save  12  cycles  of  phase  trajectory  "been  ccmplet ed? 

Arrange  for  no  counting  until  ©c  is  positive  again. 

•  • 

Is  ^  —  ©c.  negative? 

Set  up  —  ~  S>c)  . 

%  0 

Is  yOS^\  m  greater  than  preselected  £  ? 

Set  up  combined  system. 

Xs  ©c  positive? 

Is  negative? 

Is  )  positive? 

Is  0^  —  |  negative? 

Is  Oc-|  negative? 

®c.  y»©w=  01 
©e*/°©m=‘  01 

Xs  positive? 

Is  9C  ^^positive? 
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67. 


68. 


69. 


70. 


71. 


*72. 


75. 

74. 


75. 

76. 


77* 


78. 


79. 


80. 


81. 


82. 


85- 


84. 


85- 


86. 


87- 


Is  Slope  Nl  =  Ng  ? 

(l  -  Sc)  C On 


O 

9c 


-  Sf  + 

* 

Modify  FLOATEST  to  keep  system 
vhen  ec  is  negative. 

Is  0£  —  ^  ®yy\  positive? 
Store  0  L.  “  A  in  ©  yyv 


combined. 


=?  positive  value? 

Takes  off  modification 


Store  Oc  in  e  W\  . 

/* 

Is  0c  -  I  negative? 

©C 

S>C  “  OvY\  =  0? 
is  0£ 

Is  0c  -  >0  0**  positive? 

Modify  FLOATEST  to  keep  system  combined.  Takes  off  modification 

c 

when  Oc  is  positive. 

Have  ll/2  minutes  of  real  time  elapsed  since  the  last  print-out? 

Is  0c  <  1.0? 

Set  up  to  print  when  in  the  fourth  quadrant. 

Check  if  in  the  fourth  quadrant. 

Is  the  new  ©e  greater  than  the  previous 

G 

Store  the  new  time,  ©c  and  ©e. 

EXIT  (automatically  set  for  desired  jump-out),. 

Print  the  last  overshoot  stored  in  UXETBUF  (symbol  888) 

Stop  real  time  clock. 

Print  position  in  phase  trajectory  where  it  was  stopped  by  the 
real  time  clock  (symbol  999). 
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88 


.  Return  to  beginning  to  start  new  problem  if  new  parameters 
remain  to  be  evaluated.  Unconditional  stop  if  all  parameters 
have  been  evaluated. 

% 

89.  Has  the  required  number  of  overshoots  been  evaluated  prior  to 
starting  the  averaging  routine? 

90.  Add  the  <$£  just  evaluated  to  the  previous  sum  for  averaging. 

91.  Add  one  to  the  count  prior  to  averaging. 

92.  Clear  time  clock. 

o 

93.  Print  time,  and  ®c.  • 

9^.  Required  number  of  total  print-outs? 

95*  Clear  buffers. 

98.  Obtain  the  average  of  . 

97.  Subtract  1.0. 

98.  Print  the  average  of  —  1.0  and  how  many  overshoots  were 
•  utilized  to  obtain  the  average. 

99*  Is,  less  than  1.005  radians? 

100.  Clear  buffers. 

101.  Start  real  time  clock. 

102.  KX1T  (automatically  set  to  desired  jump-out). 

103*  Has  30  seconds  of  real  time  elapsed  since  the  last  print-out? 
10lf.  Print  the  position  in  the  phase  trajectory  where  it  was  stopped 
by  the  real  time  clock  and  <3:  being  less  than  1.005  radians 
(symbol  1005) • 
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Jump  Switches: 


-  For  LBUTCYC,  allows  the  maximum  overshoot  print-out  only. 

With  J  down,  a  phase  trajectory  is  printed  for  a  maximum 

« 

of  12  cycles. 

J^  -  Set  the  number  of  print-outs  depending  on  zeta.  Change 
the  number  of  overshoot  values  used  to  obtain  the  average 
value. 

J  -  Use  the  real  time  clock  for  an  automatic  recycle  to  a 

3 

new  parameter  when  have: 

a.  ll/2  minutes  maxi  mm  after  the  last  print-out. 

b.  30  seconds  maximum  and  is  1.005  radians  or  less. 


Stop  Switches:  . 

-S^  -  Stops  at  the  end  of  this  parameter  run  or  stops  at  the  end 
of  the  manual  print-out. 
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APPENDIX  B  COMPUTER  PROGRAM 


ORG 

00007 

00007 

74 

0 

00070 

40701 

fault  exf 

0 

00070 

75 

0 

SLJ 

0 

FAULPRIN 

ORG 

40000 

40000 

74 

0 

00100 

start  exf 

0 

00100 

40001 

11 

0 

1 

40004 

40623 

LAZYZETA  L^ 

0 

1 

20 

0 

40425 

STA 

0 

ZETA 

40002 

50 

1 

00000 

ENI 

1 

0 

75 

0 

40000 

SLJ 

0 

START 

40003 

75 

50 

0 

0 

SoSoo 

MANSTART  SLJ 
ENI 

0 

0 

MANUAL 

0 

40004 

12 

0 

40565 

LDA 

0 

newstop 

20 

0 

40563 

STA 

0 

STOPNOW 

40005 

12 

0 

40425 

NEWZETA  LDA 

0 

ZETA 

22 

0 

40564 

A  JP 

0 

STOP  I NDX 

40006 

32 

0 

40665 

FMU 

0 

TWO 

40007 

32 

0 

40426 

FMU 

0 

OMEGAN 

20 

0 

40404 

STA 

0 

A 

12 

0 

40624 

LDA 

0 

ZETAINDX+1 

40010 

31 

0 

40425 

FSB 

0 

ZETA 

22 

0 

40433 

A  JP 

0 

QUIKCOMP 

40011 

22 

2 

40433 

AJP 

2 

QUIKCOMP 

12 

0 

40402 

LDA 

0 

P0INT01 

40012 

20 

0 

40154 

STA 

0 

TABLES+1 

40013 

20 

20 

0 

0 

STA 

STA 

0 

0 

TABLEL+1 

TABLEM+1 

12 

0 

40661 

LDA 

0 

COUNTlO 

40014 

20 

0 

40400 

STA 

0 

tenthsec 

50 

0 

00000 

ENI 

0 

0 

40015 

75 

2 

40437 

ZZ  SLJ 

2 

BIGZETA 

50 

0 

00000 

ENI 

0 

0 

40016 

12 

0 

40652 

LDA 

0 

STOP20 

40017 

20 

12 

8 

48U2 

STA 

LDA 

0 

0 

STOPZETA 
COUNT  12 

2° 

0 

40617 

40672 

STA 

0 

COUNTZET 

40020 

12 

0 

LDA 

0 

EIGHT 

20 

0 

40620 

STA 

0 

DIVIDE 

40021 

10 

0 

00010 

ENA 

0 

10 

61 

0 

40562 

SAL 

0 

STOPRINT 

40022 

75 

4 

71000 

00000 

GOAROUND  SLJ 

4 

DECOF 

50 

0 

ENI 

0 

0 

40023 

02 

0 

40427 

02 

0 

FORZEROS 

04 

0 

00000 

04 

0 

0 

40024 

12 

0 

40664 

LDA 

0 

ONE 

20 

0 

40156 

STA 

0 

UDOT 

40025 

10 

0 

00000 

ENA 

0 

0 

20 

0 

00000 

STA 

0 

0 

40026 

20 

0 

40162 

STA 

0 

THETA 

20 

0 

40161 

STA 

0 

THETADOT 

40027 

20 

0 

40157 

STA 

0 

U 

20 

0 

40155 

STA 

0 

T 

40030 

20 

0 

40377 

STA 

0 

INDEX 

20 

0 

40615 

STA 

0 

INDEXAVE 

40031 

20 

0 

40351- 

STA 

0 

PRI NTBUF+3 

20. 

0 

40352 

STA 

0 

PRINTBUF+4 

40032 

20 

0 

40353 

STA 

0 

PRI NT8UF+5 

20 

0 

40614 

STA 

0 

LIMINDEX 

40033 

20 

0 

-40607 

STA 

0 

LIM I TBUF+2 

20 

0 

40610 

STA 

0 

LIMITBUF+3 

40034 

20 

0 

40611 

STA 

0 

LIMITBUF+4 

LIMITBUF+5 

20 

0 

40612 

STA 

0 

40035 

20 

0 

40613 

STA 

0 

LIMITBUF+6 

" 

12 

0 

40621 

LDA 

0 

PRINTSET 

40036 

20 

0 

40326 

STA 

0 

OK2PRINT+1 

20 

0 

40335 

STA 

0 

BOUNPRIN+4 

40037 

12 

0 

40622 

LDA 

0 

NEWLIMPT 

20 

0 

40544 

STA 

0 

LIMPRINT+1 

40040 

12 

0 

40253 

LDA 

0 

FLOTDATA 

20 

0 

40071 

STA 

0 

FLOATEST 

40041 

12 

0 

40254 

■  LDA 

0 

FLOTDATA+1 

20 

0 

40072 

STA 

0 

FLOATEST+1 

40042 

75 

0 

40452 

NEWINPUT  SLJ 

0 

I NDXVALU 

50 

0 

00000 

ENI 

0 

00000 

131 


40043 

40044 

40045 

40046 

40047 

40050 

40051 

40052 

40053 

40054 

40055 

40056 

40057 

40060 

40061 

40062 

40063 

40064 

40065 

40066 

40067 

40070 

40071 

40072 

40073 

40074 

40075 

40076 

40077 

40100 

40101 

40102 

40103 

*0104 

40105 

40106 

40107 


ll 

8 

48IJ8 

8 

Efhl&R 

12 

0 

40513 

LDA 

0 

NOS  I N< 

20 

il 

0 

8 

8?8Xl 

00000 

. 

STA 

fM 

0 

4 

0 

floatest 

DECOF 

0 

0 

40404 

• 

01 

0 

A 

06 

0 

00001 

06 

0 

72 

75 

0 

4 

40046 

71000 

• 

RAO 

SLJ 

0 

4 

deJof 

01 

0 

40412 

01 

0 

MOINERT 

06 

0 

00000 

06 

0 

0 

72 

0 

40510 

RAO 

0 

LINEAR+10 

75 

4 

71000 

SLJ 

4 

DECOF 

01 

0 

40420 

01 

0 

KMOTCONS 

06 

0 

00000 

06 

0 

0 

10 

0 

00000 

combined 

ENA 

0 

0 

20 

0 

40160 

STA 

0 

OX 

20 

0 

40163 

STA- 

0 

QY 

75 

•4 

60200 

SLJ 

4 

RUNGE 

00 

0 

40153 

0 

0 

TABLES 

00 

0 

40164 

0 

0 

DERIVES 

75 

4 

60201 

SLJ 

4 

RUNGE+1 

50 

0 

00000 

ENI 

0 

0 

75 

0 

40310 

SLJ 

0 

COMBPRIN 

50 

0 

00000 

ENI 

0 

0 

12 

0 

40155 

LDA 

0 

T 

20 

0 

40174 

STA 

0 

TL 

20 

0 

40210 

STA 

0 

TM 

1 2 

0 

40156 

LDA 

0 

UDOT 

20 

0 

40175 

STA 

0 

VDOT 

33 

0 

40411  ' 

FDV 

0 

RHO 

20 

0 

40211 

STA 

0 

WDOT 

12 

0 

40157 

LDA 

0 

U 

20 

0 

40176 

STA 

0 

V 

20 

0 

40200 

STA 

0 

THETADL 

33 

0 

40411 

FDV 

0 

RHO 

20 

0 

40212 

STA 

0 

W 

20 

0 

40214 

STA 

0 

THETADM 

12 

0 

40162 

LDA 

0 

THETA 

20 

0 

40201 

STA 

0 

THETAL 

33 

0 

40411 

FDV 

0 

RHO 

20 

0 

40215 

STA 

0 

THETAM 

50 

0 

00000 

ENI 

0 

0 

12 

0 

40162 

floatest 

LAC 

0 

THETA 

30 

0 

40664 

FAD 

0 

ONE 

32 

0 

40424 

FMU 

0 

OMEGANSO 

33 

0 

40161 

FDV 

0 

THETADOT 

31 

0 

40404 

FSB 

0 

A 

30 

0 

40423 

FAD 

0 

B 

22 

2 

40053 

A  JP 

2 

COMBINED 

12 

0 

40161 

LDA 

0 

THETADOT 

22 

3 

40270 

A  JP 

3 

MODIFY 

50 

0 

00000 

ENI 

0 

0 

10 

0 

00000  • 

LOAD 

ENA 

0 

0 

20 

0 

40177 

STA 

0 

OL 

20 

0 

40202 

STA 

0 

OLL 

20 

0 

40213 

STA 

0 

OM 

20 

0 

40216 

■  STA 

0 

OMM 

75 

4 

60200 

SLJ 

4 

RUNGE 

00 

0 

-40172 

0 

0 

TABLEL 

00 

0 

40203 

0 

0 

DERIVL 

75 

4 

60201 

SLJ 

4 

RUNGE+1 

50 

0 

00000 

ENI 

0 

0 

75 

4 

60200 

MOTOR 

SLJ 

4 

RUNGE 

50 

0 

00000 

ENI 

0 

0 

00 

0 

40206 

0 

0 

TABLEM 

00 

0 

40217 

0 

0 

DERIVM 

75 

4 

60201 

SLJ 

4 

RUNGE+1 

50 

0 

00000 

ENI 

0 

0 

13 

0 

40215 

CONTEST 

LAC 

0 

THETAM 

32 

0 

40411 

FMU 

0 

RHO 

30 

0 

40201 

, 

FAD 

0 

THETAL 

22 

0 

40112 

AJP 

0 

BOUNCE 

40110 

22 

3 

40112 

AJP 

3 

BOUNCE 

40111 

Si 

0 

3 

J8S19 

FSB 

AJP 

0 

3 

FL(5TPRIN 

50 

0 

00000 

ENI 

0 

0 

40112 

75 

0 

40331 

BOUNCE  SLJ 

0 

BOUNPRIN 

50 

0 

00000 

ENT 

0 

0 

40113 

12 

0 

40200 

LDA 

0 

THETADL 

20 

0 

40212 

STA 

0 

W 

40114 

12 

33 

0 

0 

40417 

40405 

LDA 

FDV 

0 

0 

RHOSO 

INERTRAT 

40115 

31 

0 

40407 

FSB 

0 

RESTITUT 

32 

0 

40200 

FMU 

0 

THETADL 

40116 

20 

0 

40200 

STA 

0 

THETADL 

12 

0 

40664 

LDA 

0 

ONE 

40117 

30 

0 

40407 

FAD 

0 

RESTITUT 

40120 

32 

0 

40411 

FMU 

0 

RHO 

32 

0 

40214 

FMU 

0 

THETADM 

40121 

30 

0 

40200 

FAD 

0 

THETADL 

33 

0 

40414 

FDV 

0 

TOTINERT 

32 

0 

40412 

FMU 

0 

MOINERT 

40122 

20 

0 

40200 

STA 

0 

THETADL 

20 

0 

40176 

x  STA 

0 

V 

40123 

13 

0 

40176 

LAC 

0 

V 

32 

0 

40423 

FMU 

0 

B 

40124 

20 

0 

40175 

STA 

0 

VDOT 

40125 

13 

0 

40214 

LAC 

FMU 

0 

THETADM 

32 

0 

40411 

0 

RHO 

30 

0 

40212 

FAD 

0 

W 

40126 

32 

0 

40407 

FMU 

0 

RESTITUT 

30 

0 

40200 

FAD 

0 

THETADL 

40127 

33 

0 

40411 

FDV 

0 

RHO 

40130 

20 

0 

40214 

STA 

0 

THETADM 

20 

0 

40212 

STA 

0 

W 

13 

0 

40212 

LAC 

0 

W 

40131 

32 

0 

40422 

FMU 

0 

D 

20 

0 

J8IH 

STA 

0 

WDOT 

40132 

13 

0 

LAC 

0 

THETAL 

30 

0 

40664 

FAD 

0 

ONE 

40133 

32 

0 

40421 

FMU 

0 

C 

30 

0 

40211 

FAD 

0 

WDOT 

40134 

20 

0 

40211 

STA 

0 

WDOT 

75 

0 

40332 

SLJ 

0 

B0UNPRIN+1 

40135 

50 

0 

00000 

COMBTEST  ENI 

0 

0 

12 

0 

40214 

LDA 

0 

THETADM 

40136 

32 

0 

40411 

FMU 

'  0 

RHO 

40137 

31 

0 

40200 

FSB 

0 

THETADL 

22 

3 

40276 

AJP 

3 

COMPLMNT 

40140 

50 

0 

00000 

ENI 

0 

00000 

65 

0 

40403 

THS 

0 

DIFFERNT 

40141 

75 

0 

40301 

SLJ 

0 

PREPLANS 

12 

0 

40200 

LDA 

0 

THETADL 

20 

0 

40161 

STA 

0 

THETADOT 

40142 

20 

0 

40157 

STA 

0 

U 

13 

0 

40157 

LAC 

0 

U  ■ 

40143 

32 

0 

40404 

FMU 

0 

A 

20 

0 

40156 

STA 

0 

UDOT 

40144 

12 

0 

40201 

LDA 

0 

THETAL 

40145 

20 

0 

40162 

STA 

0 

THETA 

13 

0 

40162 

LAC 

0 

THETA 

40146 

30 

0 

40664 

FAD 

0 

ONE 

32 

0 

40424 

FMU 

0 

OMEGANSQ 

40147 

30 

0 

40156 

FAD 

0 

UDOT 

20 

0 

40156 

STA 

0 

UDOT 

40150 

12 

0 

40174 

LDA 

0 

TL 

20 

0 

40155 

STA 

0 

T 

12 

0 

40200 

LDA 

0 

THETADL 

40151 

22 

2 

40225 

AJP 

2 

PLUSTEST 

40152 

22 

3 

40255 

AJP 

3 

MINUTEST 

75 

0 

40273 

SLJ 

0 

TEST2ER0 

50 

0 

00000 

ENI 

0 

0 

40153 

00 

0 

00000 

TABLES  OCT 

2 

00 

0 

00002 

40154 

17 

7 

15075 

'  DEC 

*01 

34 

1 

21727 

133 


40155 

88 

0 

0 

88888 

T 

DEC 

0 

40156 

20 

V 

0 

14000 

UDOT 

DEC 

1*0 

00 

0 

00000 

40157 

00 

0 

00000 

U 

DEC 

0 

00 

0 

00000 

40160 

00 

0 

00000 

QX  . 

DEC 

0 

00 

0 

00000 

40161 

88 

0 

0 

0 

ooooo 

ooooo 

ooooo 

THETADOT 

DEC 

0 

40162 

00 

THETA 

DEC 

0 

00 

0 

ooooo 

40163 

00 

0 

ooooo 

QY 

DEC 

0 

00 

0 

ooooo 

40164 

13 

0 

40157 

DERIVES 

LAC 

0 

u 

32 

0 

40404 

EMU 

0 

A 

40165 

20 

0 

40156 

STA 

0 

UDOT 

13 

0 

40162 

LAC 

0 

THETA 

40166 

30 

0 

40664 

FAD 

0 

ONE 

32 

0 

40424 

FMU 

0 

0ME6ANS0 

40167 

30 

0 

40156 

FAD 

0 

UDOT 

40170 

i? 

0 

0 

48139 

-- 

STA 

lda 

0 

0 

UDOT 

U 

20 

0 

40161 

STA 

0 

THETADOT 

40171 

75 

0 

60202 

SLJ 

0 

RUNGE+2 

50 

0 

OOOOO 

ENI 

0 

0 

40172 

00 

0 

OOOOO 

TABLEL 

OCT 

2 

00 

0 

00002 

40173 

17 

7 

15075 

DEC 

.01 

34 

1 

21727 

40174 

00 

0 

OOOOO 

TL 

DEC 

0 

* 

00 

0 

ooooo 

40175 

00 

0 

ooooo 

VDOT 

DEC 

0 

00 

0 

ooooo 

40176 

00 

0 

ooooo 

V 

DEC 

0 

00 

0 

ooooo 

40177 

00 

0 

ooooo 

QL 

DEC 

0 

00 

0 

ooooo 

40200 

00 

0 

ooooo 

THETADL 

DEC 

0 

00 

0 

ooooo 

40201 

00 

0 

ooooo 

thetal 

DEC 

0 

00 

0 

ooooo 

40202 

00 

0 

ooooo 

QLL 

DEC 

0 

00 

0 

ooooo 

40203 

13 

0 

40176 

DERIVL 

LAC 

0 

V 

32 

0 

40423 

FMU 

0 

B 

40204 

20 

0 

40175 

STA 

0 

VDOT 

40205 

12 

0 

40176 

LDA 

0 

V 

2° 

0 

40200 

STA 

0 

THETADL 

.  75 

0 

60202 

SLJ 

0 

RUNGE+2 

40206 

00 

0 

OOOOO 

TABLEM 

OCT 

2 

00 

0 

00002 

40207 

17 

7 

15075 

DEC 

*01 

34 

1 

21727 

40210 

00 

0 

OOOOO 

TM 

DEC 

0 

00 

0 

OOOOO 

40211 

00 

0 

ooooo 

WDOT 

DEC 

0 

00 

0 

OOOOO 

40212 

00 

0 

OOOOO 

W 

DEC 

0 

00 

0 

ooooo 

40213 

00 

0 

ooooo 

QM 

DEC 

0 

00 

0 

ooooo 

40214 

00 

0 

ooooo 

THETADM 

DEC 

0 

00 

0 

ooooo 

40215^, 

00 

0 

ooooo 

THETAM 

DEC 

0 

00 

0 

ooooo 

40216 

00 

0 

ooooo 

QMM 

DEC 

0 

00 

0 

ooooo 

40217 

13 

0 

40201 

DERIVM 

LAC 

0 

THETAL 

40220 

32 

0 

40421 

FMU 

0 

C 

20 

0 

40211 

STA 

0 

WDOT 

- 

13 

0 

40212 

LAC 

0 

W 
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(Numbers  indicated  in  tenths  followed  hy  power  of  10) 
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APPENDIX  D  1 

TYPICAL- PRINT  OUT  ’  ILLUSTRATING  LIMIT  CYCLE 
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2 

2$U>079999999 

lm/7^.  11111111 

fjfr  59999999 

e  99999999 

A  29999999 

P  99999999 

"3^99999999- 

1 

Tl.  89999999 

99999999 

32000000 

■fu  47999999 

p1"  99999999 

K  99999999 

K/J*99999999 

1 

fnv/Tm 32000000 

1 

fJJl  53333333 

99999999 

■j  39999999 

fuimWer  1 

-t  35400000 

1 

4t  17098208- 

1 

ec  13104190 

1 

nio>iw\ w  m  c 

11580000 

2 

17034145- 

1 

10721242 

1 

3 

23279999 

2 

18292628- 

2 

10323347 

1 

4 

38859999 

2 

35712402- 

3 

10200613 

1 

5 

58729999 

2 

72447348- 

4 

10138780 

1 

- 

6 

83939998 

2 

11243587- 

4 

10098946 

1 

7 

11619999 

3 

10685365- 

5 

10071240 

1 

8 

15810000 

3 

47854719- 

7 

10051391 

1 

9 

21077001 

3 

29068282- 

8 

10037079 

1 

1005 

215160011794  3 

-167815371770  - 

2 

999996862854 

|OaC  indicates  less  tUa*  t.005*  ^0  -sec.  o- f 


real  tv wt  \\as  elapse*}  sivsee  tVie  ta*jt 
vuas  records  Hou)  st\ocAy  V 

TV\«.  t;»v\<r  .  Velocity  t  D*d  ^  I  ci  it  cov^^t^d 

po\y>t  is 

APPENDIX  D  2 

TYPICAL  PRINT  OUT  ILLUSTRATING  NO  LIMIT  CYCLE 
t  Numbers  indicated  in  tenths  followed  by  power  of  10 ) 


genTA  7.U64  no. 19 

Backlash  in  second  order  feedback  contro 


3  2768  001  61444  9 

DUDLEY  KNOX  LIBRARY 


